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EFFECT  OF  SMALL  INITIAL  IRREGULARITIES  ON 
THE  STRESSES  IN  CYLINDRICAL  SHELLS 

1 . INTRODUCTION 

1.1.  General 

In  the  conventional  amall-defleotion  analysis  of  a thin 
circular  cylinder  under  the  action  of  uniform  normal  pressure, 
the  internal  stress  developed  in  the  tube  is  a constant  circum- 
ferential foroe  with  no  bending  moment  in  either  longitudinal  or 
circumferential  direction  at  any  point.  However,  as  is  well  Jcnowif 
(1,  2),  if  the  uniform  external  pressure  exceeds  a certain  limit 
known  as  the  critical  value,  the  cylinder  will  no  longer  maintain 
its  circular  shape.  Buckling  then  occurs.  In  the  problem  which 
is  under  investigation  in  this  manuscript  it  is  assumed  that 
the  external  pressure,  p,  is  always  less  than  this  critical  value, 
Por.>  although  the  ratio  P/Pcr,  is  not  necessarily  small  oompared 
with  unity. 

If  the  cross  seotion  of  a cylinder  deviates  from  the  cir- 
cular shape  due  to  defeots  of  manufacture , imperfeot  fabrication, 

1.  Numbers  in  parentheses,  unless  otherwise  identified,  refer  to 
the  Bibliography  at  the  end  of  this  report. 
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improper  handling,  or  the  action  of  external  forces  or  shook,  the 
strength  of  the  tube  is  considerably  reduoed  because  bending 
moment  is  introduced.  The  most  common  imperfections  in  tubes  are 
an  Initial  ellipticity  and  local  dents.  Since  the  failure  of 
cylinders  under  uniform  external  pressure  depends  very  much  upon 
the  various  kinds  of  imperfections  in  them,  it  seems  quite  impor- 
tant to  investigate  the  effeot  of  such  imperfections  on  the 
strength  of  tubes. 

When  the  initial  irregularity  is  uniform  along  the  length 
of  the  tube,  the  problem  becomes  a two-dimensional  one  and  a method 
of  analysis  for  a circular  ring  with  initial  irregularities  can  be 
applied  to  a tube  by  substituting  the  flexural  rigidity  of  a 
shell,  D = BhVl2(l  - V2)»  for  the  El  of  a ring. 

Theoretical  investigations  on  the  subject  of  a circular 
ring  with  initial  irregularities  have  been  made  by  M.  Marbeo  (3), 

S.  Timoshenko  (4),  M.  A.  Thuloup  (5)  and  G.  Salet  (6).  Id.  Marbeo 
originated  an  elegant  method  for  solving  this  problem,  using 
small-deflection  theory.  However,  when  the  uniform  external  pres- 
sure is  not  small  compared  with  its  critical  value,  the  deflections 
of  the  ring  axis  resulting  from  deformations  due  to  external  loads 
are  of  great  significance  in  the  final  result.  This  significance 
arises  from  the  fact  that  in  an  irregular  ring  the  axial  stress 
contributes  to  the  bending  moment.  This  contribution  is  ignored 
in  small-deflection  theory. 

A theory  of  thin  shells  has  been  well  established  for  many 
years.  A.  E.  H,  Love  (7)  contributed  an  especially  important  part 
of  the  literature  by  making  a mathematical  study  of  the  subject 
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and  presenting  it  in  a comprehensive  manner.  This  theory  has  been 
subjected  to  many  reviews,  one  of  the  most  recent  of  whioh  is  that 
due  to  W.  H.  Osgood  and  J.  A.  Joseph  (8).  The  expressions  for 
change  of  curvature  and  the  equations  of  equilibrium  given  in  this 
paper  differ  slightly  from  those  derived  by  Love.  However,  under 
the  assumption  of  inextensibility  of  the  middle  surface  of  the 
shell,  and  neglecting  terms  of  a non-linear  nature,  the  results 
derived  by  Love  and  by  Osgood  and  Joseph  are  identical.  Among 
other  investigators  who  have  contributed  to  the  theory  of  thin 
shells,  H.  Neuber  (9),  J.  Barthelemy  (10),  A.  L.  Goldenweiser  (11), 
V.  Z.  Vlasov  (12),  H.  L.  Langhaar  (13)  and  P.  S.  Epstein  (14)  may 
be  mentioned. 

Although  the  theoretical  background  of  the  three-dimensional 
case  has  been  developed,  the  application  of  the  theory  to  the 
problem  of  a cylindrical  shell  with  initial  imperfections  has 
scarcely  been  touched.  This  may  be  due  to  the  complicated  nature 
of  the  mathematical  expressions  for  the  relations  between  strains, 
changes  of  curvature  and  displacements,  and  equations  of  equili- 
brium. So  far  as  the  writers  know,  only  one  paper,  concerning  a 
shell  of  variable  radius,  written  by  E.  F.  Burmistrov  (15) , bears 
directly  on  the  point  in  question.  It  is  of  the  utmost  practical 
Importance,  therefore,  to  develop  a method  of  analysis,  even  if 
it  is  a very  approximate  one,  to  Investigate  the  effect  of  Initial 
irregularities  on  the  strength  of  cylindrical  shells. 

The  present  investigation  is  undertaken  with  two  objectives; 
firstly,  to  develop  a method  of  analysis,  based  on  large-defleotion 
theory,  for  a circular  ring  with  initial  out-of-roundness  and, 
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seoondly,  to  solve  a practical  three-dimensional  problem  for  a 
thin  cylindrical  shell  with  small  initial  irregularities. 

The  first  method  of  analysis  presented  in  this  report 
for  a circular  ring  with  initial  irregularities,  in  whioh  both 
large  deflection  and  extension  of  the  central  line  of  the  ring 
axis  are  taken  into  consideration,  is  a numerical  approach  whioh 
is  a combination  of  two  well-known  procedures.  A method  of  suc- 
cessive approximations  in  conjunction  with  the  numerical  procedure 
of  integration,  due  to  N.  M.  Newmark  (16)  is  used.  This  method, 
which  is  applicable  whether  or  not  the  Initial  imperfection  is 
small  and  under  any  loading  condition,  has  great  generality.  Due 
to  the  fact  that  the  procedure  involves  a considerable  amount  of 
numerical  computation,  a simplified  method  for  predicting  the 
maximum  bending  moment  developed  in  a ring  subjected  to  uniform 
pressure  is  developed.  The  simplified  method  is  approximate,  but 
it  1 s direct.  The  results  obtained  by  these  two  methods  check 
extremely  well. 

A problem  of  a cylindrical  shell  with  a dent  varying  in 
magnitude  along  the  longitudinal  axis  of  the  cylinder  is  also  solved 
in  this  dissertation  for  a uniform-pressure  load.  By  assuming 
that  the  initial  irregularity  is  small  and  that  its  square  can  be 
neglected,  all  fundamental  quantities,  such  as  curvatures,  strain 
and  stress  components,  whioh,  in  general,  have  very  complicated 
mathematical  expressions,  are  expressed  as  comparatively  simple 
forms.  Three  governing  displacement  equations  of  equilibrium 
for  the  case  of  a circular  cylindrical  shell  under  uniform  pres- 
sure are  obtained.  They  contain  correction  terms  due  to  the 
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assumed  initial  irregularities.  These  equations  are  then  solved 
and  the  predictions  of  the  analysis  investigated  from  an  engineer- 
ing point  of  view. 

1.2.  Acknowledgment 
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Department  of  the  Navy,  under  contract  N6ori-71,  Task  Order  VI. 
This  report  is  based  on  a dissertation  submitted  to  the  Graduate 
Faculty  of  the  University  of  Illinois  by  T.  S.  Wu,  in  partial 
fulfillment  of  the  requirements  for  the  degree  of  Doctor  of 
Philosophy  in  Engineering,  1952. 
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2.  CIRCULAR  RING  WITH  SMALL  INITIAL  IRREGULARITIES 

2.1.  Recapitulation  of  Small-Def leotlon  Theory 
2.1.1.  Method  of  Analysis 

In  this  theory  several  assumptions  are  made: 

(1) .  The  cross-sectional  dimensions  of  the  ring  are  assumed 
to  be  small  in  comparison  with  the  unstrained  radius  of  curvature. 
The  stress  distribution  is,  on  any  cross  section  then,  nearly  a 
linear  one. 

(2) .  The  central  line  of  the  ring  is  regarded  as  inextensible. 
Any  extension  of  the  central  line  which  might  be  produced  by  the 
given  forces  would  be  extremely  small  compered  with  the  flexure 
actually  produced. 

(3) .  The  uniform  external  pressure  must  be  much  less  than 
the  critical  value  and,  therefore,  the  effect  of  the  axial  com- 
pressive force  on  the  bending  can  be  neglected. 

(4) *  The  deformations  and  displacements  due  to  the  action  of 
external  loads  are  very  small  oompered  with  the  dimensions  cf  the 
ring  concerned.  Then  the  same  dimensions  can  be  used  in  the  un- 
strained as  in  the  strained  state  of  the  ring  and  the  total 
deflections  resulting  from  suoh  deformations  will  usually  be  tos 
small  to  modify  the  shape  of  the  ring  materially. 

(5) *  The  residual  or  initial  stress  due  to  the  source  of  the 
imperfection  is  entirely  disregarded. 

Marbeo’s  method,  introduced  in  1908,  is  the  simplest  and 
the  most  convenient  one  for  the  solution  of  the  problem  under 
investigation.  It  has  also  been  derived,  by  R.  Mayer  (17)  and 
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H.  Neuber  (18),  using  the  method  of  aroh  analysis.  It  leads  to 
the  oonolusion  that  in  the  oase  of  a ring  of  any  shape  and  of 
variable  thickness,  under  the  aotion  of  uniform  pressure,  the 
internal  resisting  forces  and  the  internal  bending  moment  acting 
at  any  point  (x,  y)  on  the  ring  are  given  by  the  expressions 
(Fig.  2. 1.1.1) s 

H = -p(y  - b) , V = p(x  - a)  2. 1.1.1 

and  lfi  - p(rc2  - fc*  )/2  2. 1.1. 2 

where  £c2  s(x  - a}2  -f  (y  - b)2 

The  symbols  a and  b represent  the  coordinates  of  the  cen- 
ter of  a "node"  circle  of  radius  rc  , and  is  the  distance  from 
the  point  x,  y to  the  point  a,  b.  All  coordinates  are  measured 
with  respect  to  principal  axes  X,  Y with  origin  at  the  centroid. 

It  can  be  seen  that  the  bending  moment  at  any  point  on 
the  ring  lying  outside  the  node  circle  is  positive;  for  points 
inside  the  node  circle,  bending  moment  is  negative.  The  bending 
moment  at  any  point  on  the  node  circle,  i.e.  = T*c  , is 
equal  to  zero;  so  that  the  node  circle  represents  the  funicular 
curve  of  the  ring.  The  resultant  of  H and  V is  equal  to  px^c 
and  acts  in  a direction  perpendicular  to  • 

The  three  unknown  quantities  a,  b and  rc  oan  be  easily 
determined  if  the  node  oircle  is  located.  The  coordinates  of 
the  center  of  the  node  oircle  with  reference  to  the  principal  axes 
of  a given  shape  of  ring  are: 

a = [J(x2  + y2)xds/El]/2[Jx2ds/El] 
b~[J(x2  f y2)yds/El]/2[Jy2ds/El] 


2.1.1. 3 
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The  radius  of  the  node  oirole  la  given  by: 

r02  = (.2  ♦ b2)  + £/(*2  + y2)d./Bl]  /[Jd./Sl]  2.1J.4 

All  integrals  are  to  be  taken  over  the  entire  circumference  of  the 
ring. 

It  is  evident  that  if  the  shape  of  a ring  is  symmetrical 
with  respect  to  the  z or  y axis,  a or  b,  respectively,  will  be 
equal  to  zero  and,  if  the  ring  has  two  axes  of  symmetry,  the  oen- 
ter  of  the  node  circle  C will  coincide  with  the  centroid  G of  the 
ring. 

When  the  shape  of  a ring  is  known,  the  centroid  of  the 
given  shape  and  its  principal  axes  with  reference  to  any  arbitrary 
axes  X'  and  Y'  can  be  determined  by  the  well-known  expressions: 

* a - [J*'d./Ei]/(Jd«/EiJ 

AeieXey 

7q  - [fy'da/El]/[Jd»/El] 

and  tan2V  » [2]l'y'da/?.l]/^[j(l'  )!<l8/El]-[j(y'  )2d*/ElJj 

where  * angle  between  X and  X7  axes. 

From  the  known  coordinates  (x,  y)  of  any  point  on  the  ring 
with  reference  to  its  principal  axes  X and  Y,  the  three  quantities 
a,  b,  and  r0  in  Equations  (2. 1.1. 3)  and  (2. 1.1. 4)  oan  be  evaluated 
by  performing  the  integrations  either  analytically  or  numerically. 
Beoause  of  the  complicated  given  shape  of  a ring,  difficulty  may 
be  encountered  in  performing  the  integration  of  these  quantities 
analytically.  For  this  reason  the  numerical  integration  procedure 
due  to  N,  M.  Newmark  (16),  whioh  yields  quite  aoourate  results 
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with  a reasonable  amount  of  work,  is  to  be  preferred  in  the  major- 
ity of  oases. 

If  the  initial  irregularity  is  very  small,  a further  sim- 
plification oan  be  made. 

M p [ (rc  + fc  )/2]  (re  - fc  ) 

Since  the  initial  irregularity  is  very  small,  rc  and  fc  differ 
from  R0  only  by  a very  small  amount  and,  therefore,  it  oan  be 
assumed  that  (rc  + fc  )/2  is  approximately  equal  to  R 0 • Then, 

M pR 4 ( r ) 2. 1.1. 6 

Slnoe  the  axial  foroe  in  a perfect  circular  thin  ring  is  pR„, 
Equation  (2. 1.1. 6)  indicates  that  for  small  initial  irregularities 
the  bending  moment  at  any  point  is  equal  to  this  foroe  multiplied 
by  the  irregularity  at  thst  point  measured  with  respect  to  the 
node  circle.  The  bending  moment,  then,  bears  a linear  relation- 
ship to  the  irregularity  measured  with  reference  to  the  node  oir- 
cle.  The  area  between  the  node  circle  and  the  ring  axis  represents 
the  bending  moment  diagram  for  the  ring.  From  a geometric  point 
of  view  the  area  bounded  by  the  ring  axis  and  the  funicular  curve 
(which  is  a cirole)  balanoe  eaoh  other* 

2.1.2.  Illustrative  Problems 

To  illustrate  the  application  of  the  Marbeo  method  to  a cir- 
cular ring  with  initial  Irregularities  under  the  action  of  a uni- 
form external  pressure,  two  numerical  problems  are  worked  out.  The 
first  is  the  oase  of  an  elliptical  ring  while  the  second  concerns 
a olroular  ring  having  two  symmetrical  looal  dents. 


(1).  Elliptical  Ring 

The  equation  of  an  elliptical  ring  is  expressed  by: 
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x2/(.99R0  )2  ♦ y2/{1.01R  * )2  ~ 1 2. 1.2.1 

Sinoe  the  ellipse  has  two  axes  of  symmetry,  the  center  of  the 
node  circle  and  the  centroid  of  the  ellipse  coinoide  (Fig.  2. 1.2.1). 
Both  a and  b are  zero  and  Equation  (2. 1.1. 4)  beoomes 

rc2  = [J(x2  + y2)ds/EI  ] /(Jds/BlJ  2. 1.2.2 

In  order  to  simplify  the  numerical  integration,  Equation 
(2. 1.2.1)  is  transformed  from  Cartesian  into  polar  coordinates. 

f2,  = (l.OIRo  )2/(k2cos2<p  + sin2^  ) 2. 1.2.3 

rc2  = [ / B ^d^/El]  /[fBdf/El]  2. 1.2. 4 

k a*  1.01/. 99 

B = [ 

Because  of  the  double  symmetry  only  one  quarter  of  the 
ring  need  be  considered.  To  illustrate  the  numerical  method  six 
equal  angular  divisions  are  taken  in  the  quarter  ring.  It  is  of 
advantage  to  take  all  the  divisions  equal  when  possible,  however, 
this  1 8 not  necessary.  The  numerical  values  of  p and  r^  com- 
puted  from  Equation  (2. 1.2. 3)  and  (2. 1,2. 4)  respectively  are 
given  in  Table  (2. 1.2.1).  The  maximum  axial  force,  bending  moment 
and  the  fiber  stresses  due  to  this  maxium  axial  force  and  bend- 
ing moment  are  given  in  Table  (2. 1.2.4) • 

{2. ) Circular  Ring  with  Two  Symmetrical  Local  Dents 

A circular  ring  with  two  symmetrical  local  dents  is  shown 
in  Fig.  (2. 1.2. 2).  The  dented  portion  for  -<j < <p0 


and 

in  which 
and 


and 
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JH6<  + is  defined  by 

^ = (1  - & cos3  t T<f  /2f0  ) 2. 1.2.5 

and  for  other  values  of  (f  , ^ • R*  . 

Here  S and  (f0  are  two  parameters  which  determine  the  depth  of 
the  dent  and  the  angular  area  over  which  it  extends. 

As  before,  the  ring  has  two  axes  of  symmetry,  a and  b are 

equal  to  zero  and  r 2 is  given  by  Equation  (2. 1.2. 4}.  One  quar- 

c 

ter  of  a ring  with  six  equal  angular  divisions  is  again  used.  In 
order  to  take  the  effect  of  the  dent  into  consideration,  ten  sub- 
divisions are  provided  in  the  dented  portion.  For  two  oases 

3 15'  , s - 0,01  15"*  S =■  0.1  th«  numerical  values 

of  ^ and  r0  computed  from  Equation  (2. 1.2. 5)  and  (2. 1.2. A)  are 
given  in  Tables  (2. 1.2. 2)  and  (2. 1.2. 3)  respectively.  The  maxi- 
mum axial  force  and  bending  moment  and  the  fiber  stresses  for 
these  two  cases  due  to  this  maximum  axial  force  and  bending 
moment,  are  given  in  Table  (2. 1.2. 4). 

Examining  the  numerical  results  of  these  two  illustrative 
problems  shown  in  Table  (2. 1.2. 4)  it  is  interesting  to  note  that, 
although  the  initial  irregularity  of  the  ring  is  very  small  (1/100 
of  the  radius),  the  maximum  fiber  stress  due  to  bending  moment 
is  from  three  to  six  times  as  great  as  that  due  to  axial  force. 
The  maximum  total  fiber  stress  in  a ring  with  an  initial  irregu- 
larity of  this  magnitude  is  magnified  four  to  seven  times  that  in 
a perfect  circular  ring  for  R./h*50  or  100.  If  the  initial 
irregularity  increases  to  1/10  of  the  radius,  the  maximum  total 
fiber  stress  is  from  28  to  54  times  that  in  a perfect  circular 
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ring.  This  shows  the  great  importance  of  Initial  irregularity 
from  the  point  of  view  of  safety  in  design. 

The  formula  for  bending  moment  given  by  Equation  (2. 1.1. 2) 
is  exact.  The  approximate  formula  for  bending  moment  given  by 
Equation  (2. 1.1. 6)  is  only  applicable  to  cases  in  which  the  ini- 
tial irregularity  is  very  small.  When  the  initial  irregularity 
of  a ring  ranges  from  1/100  to  1/10  of  the  radius,  the  dis- 
crepancy between  the  * exact  * bending  moment  and  the  approximate 
one  varies  from  0.5%  to  5%.  For  most  practical  cases,  therefore, 
the  approximate  formula  (Equation  2. 1.1. 6)  can  be  used  without 
introducing  serious  error. 

In  the  above  mentioned  illustrative  problems,  it  has  been 
assumed  that  the  ring  under  investigation  has  a uniform  cross 
section  so  that  the  quantity  El  is  a constant  and  oan  be  taken  as 
unity  in  numerical  calculations.  If  the  thiokness  of  a ring 
varies  along  its  circumference,  the  average  value  of  El  for  every 
division  in  the  numerical  integration  should  be  evaluated.  The 
values  of  B and  B ^ divided  by  the  corresponding  values  of  El 
for  the  different  divisions  are  then  to  be  used  in  the  integration. 
All  other  steps  remain  unchanged. 
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re/B. 

1.000000 

0.999296 

0.992901 

M/pR* 

0.010050 

0.009246 

0.087926 

Approx.  M/pR* 

0.010000 

0.009296 

0.092900 

B 

50 

150.8 

138.7 

1319 

P 

100 

603.0 

554.8 

5276 

B 

Jk 

50 

50.50 

50.00 

50 

P 

h 

100 

101.0 

100.0 

100 

Z<r 

50 

4.0 

3.8 

27.3 

100 

7.0 

6.5 

53.7 

(r  = Tiber  street  due  to  mudmm  bending  aoaent 
<TM  = Tiber  street  due  to  aaxiana  axial  force 
2 (T  = Total  aaxlnaa  fiber  street 

pR0/h  * Tiber  stress  due  to  axial  force  in  perfect  circular  rings 


TABU  2. 1.2.4  MAXIMUM  STBXSSSS  II  BUM 
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2.2.  Large-Deflection  Theory 
2.2.1.  Method  of  Analysis 

The  method  of  analysis  presented  in  this  section  for  the 
stress  analysis  of  a circular  ring  with  initial  irregularities  sub- 
jected to  uniform  normal  pressure  is  a numerical  approach  by  means 
of  which  both  large  deflections  and  extension  of  the  ring  axis  can 
be  taken  into  consideration.  It  is  also  applicable  to  a ring  of 
variable  thickness  under  any  loading  condition. 

In  this  investigation  it  is  not  assumed  that  the  uniform 
normal  pressure  is  small  compared  with  its  critical  value.  Two 
assumptions  (1)  and  (5)  of  Article  (2.1.1)  are  still  used,  how- 
ever, in  the  present  disoussion. 

When  a ring  is  acted  on  by  a system  of  external  forces, 
there  will  always  be  a certain  displacement  of  the  ring  axis. 

The  relation  between  loads  and  displacements  can  be  made  clear  by 
consideration  of  the  physical  behavior  of  the  ring  axis  during 
deformation.  The  external  forces  induce  in  the  ring  internal 
stresses  which,  in  turn,  cause  rotations  as  well  as  axial  and 
shearing  deformations  throughout  the  ring  axis.  A deformation  in 
any  part  of  the  ring  axis  is  accompanied  by  a corresponding  dis- 
placement of  the  ring  axis,  these  two  being  interconnected  by  geo- 
metric relations.  The  effect  of  shear  on  displacement,  in  the 
case  of  a thin  ring,  is  very  small  and  will,  therefore,  be 
neglected  in  all  the  subsequent  discussions.  The  d is placements  of 
the  ring  axis  modify  its  original  shape  and  consequently  change 
the  line  of  application  of  external  loads.  In  a stable  situation, 
the  ring  will  reach  a final  equilibrium  state. 
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In  large-defleotion  theory,  the  deflection  of  the  ring  axis 
due  to  lengthening  or  shortening  must  be  considered  in  setting  up 
equations  of  equilibrium  and  equations  of  geometric  continuity. 
Before  the  ring  can  be  analyzed  according  to  large-defleotion 
theory,  it  is  therefore  necessary  to  set  up  the  equations  for 
deflections  in  terms  of  the  known  and  unknown  forces  causing  the 
deflection. 

( 1) . Relation  Between  Deformation  and  Stresses 

Fig.  (2. 2. 1.1)  shows  an  element  of  a planar  curved  ring 
in  its  final  or  deformed  state.  Under  the  action  of  external 
loads  Internal  stresses,  which  are  developed  in  the  fibers,  can 
be  represented  by  a bending  moment  M,  a normal  foroe  N acting 
through  the  centroldal  axis  on  a cross  section  normal  thereto  and 
a shearing  foroe  Q,  acting  parallel  to  the  cross  section. 

The  deformations,  which  are  produced  by  the  bending  moment 
M and  normal  foroe  N,  consist  of  an  extension  of  the  ring  axis 
and  a rotation  of  the  cross  section.  The  bending  moment  is  taken 
as  positive  when  it  decreases  the  initial  curvature  of  a ring  and 
the  normal  force  is  positive  when  it  is  tension.  The  angle  change 
and  deformation  of  the  ring  axis  are  considered  as  positive  in 
the  same  sense  as  positive  moment  and  normal  force  respectively. 

In  the  case  of  a thin  ring,  the  stress  distribution  over 
the  cross  section  of  the  ring  is  nearly  a linear  one.  Under 
these  conditions,  the  centroldal  axis  coincides  with  the  neutral 
axis  of  the  ring.  The  extension  of  the  ring  axis  due  to  rotation 
of  the  cross  section  produced  by  the  bending  moment  is  zero.  The 
extension  of  the  ring  axis  due  to  normal  foroe  is  given  by 
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A da  » 6 ds  = Nds/AE  . 2. 2. 1.1 

The  rotation  of  the  oross  section  due  to  the  extension  of 
the  ring  axis  oan  be  expressed  by 

( A d0  ^ = - € ds/R*  =■  -Nds/AER1  . 

The  rotation  of  the  cross  section  due  to  bending  moment  is 
( A d0  )2  =s  lids'  /El  . 

The  total  rotation  is , then,  equal  to 

Ad0  - (lids'/EI)  - ( 6 ds/R' ) . 2. 2. 1.2 

The  curved  length  ds'  and  final  curvature  1/R'  of  an  element  in 
the  deformed  state  can  be  found  by  following  relations: 

ds'  = (1  + € ids 

and  ( 1/R* ) - (l/R)  = (14/BI)  - (t/R'i  2. 2. 1.3 

Due  to  the  fact  that  the  value  of  6 is  very  small  com- 
pared with  unity,  an  approximation  which  will  not  introduce  appre- 
ciable error  can  be  made  by  neglecting  terms  of  higher  order  in  £ . 
Substituting  Bquations  (2. 2. 1.3)  into  (2. 2. 1.2)  we  have 

Ad e = (Mds/BI)  - ( 6 ds/R) . 2.2. 1.4 

(2).  Relation  Between  Displacements  and  Deformations 
In  Fig.  (2.2.1.2a),  JJ®  represents  the  initial  or  unde- 
formed state  of  an  element  of  the  ring  axis.  When  it  is  acted  by 
external  loads,  JJ#  will  be  deformed  and  displaced  to  J'J#’. 

Due  to  the  change  in  length  Ads  and  the  change  in  slope  AB  of 
the  ring  axis  caused  by  deformation,  there  will  be  corresponding 
changes  d^  and  dv)  in  the  horizontal  and  vertioal  projection  of 
the  element  respectively.  As  the  deformations  are  small  quantities 
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their  effeots  on  the  total  displacement  may  be  considered  indepen- 
dently and  they  can,  therefore,  be  separately  treated. 

Let  element  JJ o be  shifted  toward  J'J0'  until  points  J 
and  J'  coincide  with  each  other  (Fig.  2,2.1.2b).  Let  d^,  and 
dij,  be  the  deflections  of  element  JJ«  due  to  extension  of  ring 
axis.  Then, 


cos0Ads  — fccoeO  ds 


2. 2. 1.5 


dTJi  = sinG  Ads  « 0sin0  ds 

Let  d Ijz  and  dvj&  be  the  deflections  of  J J©  due  to  rotation  of 

the  element.  When  is  small,  it  can  be  assumed  that 
is  similar  to  . Then, 


dja.  = - A0(dy  + d^  ) 
dTJz  » A0  (dx  f dj  ) 

The  total  defleotions,  then,  are  given  by 


2. 2. 1.6 


d$  =*  - A0  (dy  i dty  ) + € oos0  ds 
d V}  = A0  (dx  + d^  ) + £ sinQ  ds  . 


2. 2. 1.7 


Consequently  the  rates  of  change  of  total  deflection  of  the  element 
JJ©  along  ds  are: 


d$  /ds  = - A0  (sin©  + dij/ds)  f ( cosd 

2. 2. 1.8 

d*)  /ds  * A0  (cos©  ♦ d^/ds)  4*  6 sin0 
Solving  the  above  two  equations,  we  have 

d^  /ds  » £ -(1  + * ) A0  sin0  - [(  A0  )2  - cose}/[l  f ( Ad  )*J 

dlj/ds  - { (1  + € ) A0cos0  - [(  A0  )2  -€]  sin0}/[l  f (A8  )*]. 


— 2.2. 1.9 
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The  above  two  equations  furnish  the  necessary  relations  between 
displacei&ents  and  strains  and  rotations  of  an  element  of  the  ring 
axis.  The  term  ( A 0 )£  which  is  usually  disregarded  in  the 
small- deflect! on  theory,  will  be  included  in  the  present  Investiga- 
tion. In  inextenslble  theory,  t is  taken  equal  to  zero. 

(3).  General  Governing  Equations 

When  a ring  is  subjected  to  external  loads,  it  will  deform 
and  displace  to  its  strained  state.  It  is  obvious  that  all  exteov 
nal  loads  and  internal  forces,  acting  on  any  small  element,  must 
be  in  equilibrium  and,  furthermore,  the  displacements  at  any 
point  due  to  the  deformations  must  fulfill  conditions  of  geome- 
tric continuity  (conditions  of  compatibility). 

The  equations  of  equilibrium  can  be  obtained  by  statics,  but 
in  order  to  establish  the  equations  of  compatibility  the  elastic 
deformation  of  the  ring  must  be  considered. 

The  rotation  and  displacements  at  a point  J0'  measured  with 
respect  to  a reference  point  J(  * are  given  by  (Fig.  2, 2. 1.3): 

ASj;  = £>(  ) 

3*'  ~ [-A0  (dy  f di)  ) + ecos6  dsj  2.2.1  .K) 

yjj,'  =•/**  [A0  (dx  + d£  ) + 6 sin0  ds] 

On  Integrating  the  last  two  of  Equations  (2.2.1.10)  by  parts, 
they  simplify  to: 

5*;  3 + Tj  ) " (y*  * *)o  )]d(A0  ) + Jj**  cos 6 ds 

- (x  + 5 )]d(Afl  ) + J*<sinfl  ds 

2.2.1.11 
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Since  the  ring  is  a closed  one,  points  ' and  J0*  can  be 
chosen  so  as  to  be  the  same  point.  In  order  to  provide  geometri- 
cal continuity,  the  rotation  and  displacements  at  J»'  must  be 
equal  to  zero.  Substituting  Equations  (2. 2. 1.2)  into  the  first 
equation  of  (2. 2. 1.9)  and  (2.2.1.11)  and  equating  the  rotation 
and  displacements  to  zero,  we  have 

A®-  Jm-H-  - J 

% » JM  [ty  +>))  - j + Jec«secis  * o 

2.2.1.12 


The  above  three  equations  represent  the  conditions  of  compatibil- 
ity. 

Let  the  ring  now  be  cut  at  point  J0  ’ . Internal  foroes  H., 
W,  and  Mf(  must  exist  at  J®  ’ to  produce  continuity  at  that  point 
(Fig.  2.2.1. 3).  The  bending  moment  at  any  point  J’  can  be 
written  as 

M-  Mo + t- AW1»  . 2.2.1.13 


In  the  case  of  a ring  subjected  to  uniform  external  pressure,  the 
statio  moment  and  its  increments  due  to  deflection  of  the  ring 
axis  will  be: 


yns* 


2. 2.1.14 
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By  substituting  Equations  (2.2.1.13)  into  (2.2.1.12)  they  become 
+ J"'>  El  +JiWs  Ef-  1 1C~  0 

% - + M./fty-  y.Mr’Wfff 

+J«»[(y-jy+(,f-,h)}^  -/‘[tj'yoHri"3T+^c“8d5  =0 

Vl|[(x.-X}t(  J„~  I j]  If  + Jws  [(*'-*)  H ?.•  f yj  Ej- 

2.2.1.15 

Equations  (2.2.1,15),  which  furnish  the  conditions  of 
geometrical  continuity  at  any  point,  contain  three  unknown  forces, 
H#,  V«  amd  M#  , and  two  unknown  displacements,  $ and  . 
Equations  (2.2.1. 9)  furnish  additional  relations  between  forces 
and  displacement  components.  These  five  Independent  equations 
are  sufficient  to  determine  the  five  unknowns.  Any  forces  and 
deflection  configurations  which  satisfy  Equations  (2.2.1.15)  and 
(2. 2. 1.9)  will  automatically  satisfy  the  equations  of  equili- 
brium  and  compatibility  equations.  It  can  be  seen  that  if  the 
displacements  \ and  -tj  and  the  extension  £ of  the  ring  axis 
are  neglected  in  Equations  (2.2.1.15)  these  equations  reduce  to 
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well-known  aroh  formal as. 

Equations  (2.2.1.15)  and  (2. 2. 1.9)  which  take  the  extent 
sion  and  defleotion  of  the  ring  axis  into  consideration  are 
applicable  to  a ring  of  any  shape  and  subjected  to  any  loading 
condition.  In  the  present  investigation,  only  uniform  normal 
pressure  is  of  interest. 

When  the  uniform  pressure  acting  on  the  ring  is  internal, 
a negative  value  of  p is  to  be  used  in  Equations  (2.2.1.15). 

From  equations  (2.2.1.15),  it  can  be  seen  that  all  the 
Integrals  have  physical  meanings.  In  the  first  equation, 

J ds/EI  represents  the  angle- change  between  two 
points  due  to  M,  = 1 

J [(X.-*  )+(f.-*>]ds/BI  represents  the  angle -change  between 

two  points  due  to  Ve  = 1 

j[CV_y*) -H  represents  the  angle-change  between 

two  points  due  to  H,  = 1 

j ) ds/EI  represents  the  angle-change  between 

two  points  due  to  m$ 

Jfc  ds/R  represents  the  angle-change  between  two 

points  due  to  €.  (angle-change  5 change  in  slope) 

The  quantities  in  the  second  and  third  equations  represent  the 
horizontal  and  vertical  components  of  defleotion  between  two 
points  due  to  forces,  static  moment  and  the  extension  of  the  ring 
axis. 


Due  to  the  complexity  of  these  two  sets  of  equations,  it 
is  difficult  to  obtain  an  analytical  solution  for  a general  case. 
For  this  reason,  a numerical  approach  which  serves  as  a useful 
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tool  to  attack  problems  of  this  nature  is  suggested. 

(4).  Numerical  Method 

The  method  presented  in  this  dissertation  for  the  stress 
analysis  of  a circular  ring  with  initial  out-of-roundness  under 
uniform  normal  pressure  analyzed  by  large-deflection  theory  is 
a numerical  procedure  which,  in  turn,  is  a combination  of  two 
well-known  procedures,  namely,  the  numerical  procedure  of  inte- 
gration and  the  method  of  successive  approximations  due  to  N.  M. 
Newmark  (16). 

The  general  processes  of  analysis  can  be  summarized  in 
the  following  steps: 

(i) .  Divide  the  ring  into  a number  of  equal  divisions 
along  the  curved  length.  As  mentioned  previously,  in  order  to 
simplify  the  numerical  integration,  it  is  preferable  to  have 
equal  angular  divisions. 

(ii) .  Svaluate  all  physical  quantities  such  as  slope, 
curvature,  curved  length  factor  ( ds/d ),  static  moment,  etc. 
at  all  division  points.  They  will  be  encounted  in  the  numerical 
integration.  Most  of  these  quantities  can  be  readily  calculated 
from  the  given  equation  of  the  ring.  The  static  moment  depends 
upon  the  loading  condition  as  well  as  the  geometrical  shape  of 
the  ring. 

(iii) .  Assume  two  deflection  configurations,  $ and  yj  , 
and  negleot  € .in  the  first  approximation. 

(iv) .  All  integrals  in  Equations  (2.2.1.15)  can  now  be 


evaluated  numerically. 
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(v) .  Solve  for  forces  H0,  V0  and  bending  moment  M 0 by 
means  of  Equations  (2.2.1.15)* 

(vi) .  From  H#  and  V#  , compute  N and,  in  turn,  € at  all 
division  points.  Evaluate  J*€ds/R. 

(vii) .  Calculate  the  total  angle- change , a6  , by  integrat- 
ing Equations  (2. 2.1.9)*  In  solving  actual  problems  the  total 
angle-change  at  any  two  specified  points  can  be  obtained  by  sum- 
ming up  all  components  of  angle  change,  due  to  forces,  static 
moment  and  the  extension  of  the  ring  axis,  from  the  Integrals  in 
the  first  equation  of  (2.2.1.15)  which  have  been  evaluated  in 
(iv) . 

(vii).  Two  derived  deflection  configurations,  ^ and  *j  , 
can  be  obtained  by  integration  from  Equation  (2. 2. 1.9). 

(ix).  Repeat  the  procedures  until  the  derived  and  assumed 
deflection  configurations  agree  with  each  other. 

In  the  case  of  a ring  symmetrical  about  one  axis  (which 
is  usually  taken  as  one  of  the  coordinate  axes)  only  half  of  the 
ring  has  to  be  considered.  The  horizontal  shearing  forces  at 
the  two  ends  of  the  axis  of  symmetry  are  equal  to  zero.  Since 
the  vertical  tangents  at  two  ends  remain  vertical  after  deforma- 
tion, the  relative  angle- change  at  one  end  with  respect  to  the 
other  is  zero.  Furthermore,  the  relative  vertical  displacement 
of  one  end  with  respect  to  the  other  is  also  zero.  Two  equations 
of  geometrical  continuity  are  sufficient  to  solve  for  the  two 
unknowns. 

When  the  ring  has  two  axes  of  symmetry  only  one  quadrant 
need  be  taken  into  consideration.  Due  to  symmetry  the  shear  at 
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both  ends  must  be  zero  and  the  axial  force  can  be  expressed  in 
terms  of  the  known  ordinate  and  unknown  deflection  of  either  end. 
The  only  unknown  remaining  is  the  bending  moment  whioh  can  be 
determined  from  the  condition  that  the  angle- change  at  one  end 
with  respect  to  the  other  end  is  zero. 

The  accuracy  of  the  numerical  integration,  which  deter- 
mines the  accuracy  of  the  final  result,  depends  greatly  upon  the 
nature  of  the  function  to  be  integrated  and  the  number  of  divi- 
sions chosen  along  the  ring  axis.  The  more  divisions  chosen,  the 
more  accurate  will  the  final  result  be.  On  the  other  hand  more 
numerical  computations  are  required.  For  most  practical  problems, 
six  divisions  in  one  quadrant  of  a ring  will  give  reasonably 
accurate  results.  The  equivalent  concentration  formulae  due  to 
N.  M.  Newmark  (16)  will  yield  quite  accurate  results  with  a 
reasonable  amount  of  work. 

The  rate  of  convergence  of  the  method  of  solution,  which 
is  an  iterative  procedure,  depends  on  the  correctness  of  the 
assumed  deflection  configuration  and  on  the  magnitude  of  the  exter- 
nal pressure.  A deflection  configuration  which  can  be  used  to 
hasten  the  convergence  of  the  solution  is  obtained  by  multiplying 
the  deflections  determined  by  the  small-deflection  theory  by  the 
magnification  factor  1/(1  - p/ptr  ).  The  same  technique  can  be 
applied  to  a ring  under  internal  pressure  in  which  case  the  fac- 
tor is  1/(1  + p/Pcr  )•  For  301,10  problems,  the  assumed  deflec- 
tions obtained  in  this  manner  will  be  very  close  to  the  correct 
ones.  If  the  external  pressure  is  much  less  than  its  critical 
value,  the  solution  will  converge  rapidly  whereas,  if  the  external 
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pressure  approaches  it  critical  value,  the  rate  of  convergence 
is  slower,  as  would  be  expected.  For  most  cases  a very  accurate 
result  can  be  achieved  through  four  oycles  of  iteration. 

2.2.2.  Illustrative  Problems 

In  order  to  illustrate  the  application  of  this  numerical 
method  to  practical  problems,  circular  rings  with  various  kinds 
of  imperfeotion  are  analyzed  in  this  article. 

In  large-defleotion  theory,  the  magnitude  of  the  uniform 
pressure  acting  on  the  ring  must  be  specified.  We  consider  the 
ring  to  be  subjected  to  the  action  of  a uniform  external  pres- 
sure of  magnitude  p = 1.5El/R*  for  all  illustrative  examples. 

This  implies  P/Pcr>~°«5. 

(1)  Elliptical  Ring 

The  elliptical  ring,  with  ellipticity,  e,  equal  to  1 • 
(0.99/1. 01), illustrated  in  Article  (2.1.2)  (Fig.  2. 1.2.1)  will 
be  more  accurately  analyzed  using  large-deflection  theory. 

Since  the  ring  has  two  axes  of  symmetry,  only  one  quad- 
rant, taken  with  six  equal  angular  divisions,  need  be  investi- 
gated. It  is  evident  that  the  tangents  at  division  points  0 and 
6 will  remain  vertical  and  horizontal  respectively  after  strain- 
ing. The  total  angle-change  at  point  0 with  respect  to  point  6 
is  therefore  equal  to  zero.  Then, 


M.J& + V.j'[Oc.-^*c5.-5)]#+/m.Tr+/an1* eT  - J# 


O 
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From  statics,  the  following  relations  can  be  written: 

Ho=  0,  H6  =*  -p(y6  ) 

and  VD  =»  p(x0  + 

The  integrals  J ds/EI,  Jmsds/EI  and  J (x0  - x)ds/EI  in 
Equation  (2.2.1.16)  are  constants  and  the  others  are  function  of 
^ and  yj  . Once  the  deflections  ^ and  are  either 
assumed  or  known,  the  forces  H^,  VQ  and  bending  moment  Mq  can  be 
easily  found. 

The  slope  of  the  tangent  at  any  point  on  the  ellipse  (Fig. 
2. 2. 1.4)  is  given  by 

tan©  » dy/dx  = -(O.99/1.01)2x/y  . 

The  values  of  sinC  and  cos8  at  division  points  whioh  are  required 
in  integrating  Equation  (2. 2. 1.9)  can  be  evaluated  readily. 

In  case  the  extension  of  the  ring  axis  is  to  be  included 
in  the  investigation,  one  more  numerical  constant  must  be  com- 
puted.  The  computed  constant  used  in  the  present  example, 
pR./EA  =s  0.0005 , is  based  on  these  data:  p sx  30  lbs  per  in., 
Ro=50  in.,  h =»  1 in.  and  E 30  x 10^  lbs  per  sq.  in.  In  the 
case  of  a thin  ring,  the  stability  effect  usually  predominates 
and  the  effect  of  extension  of  the  ring  axis  is  extremely  small 
in  comparison  to  that  produced  by  flexure. 

The  initial  curvature  at  any  point  is  given  by  the 
expression 


t'/R ) = [ f* zcp1-  f 0]/[ 


The  values  of  1/R  at  division  points  are  necessary  for  performing 
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the  integration  j£ds/R.  The  axial  foroe  N for  use  in  comput- 
ing £ (Pig.  2. 2. 1.4)  is  given  by  N=  -Vain©  + Hcos0  . 

As  previously  mentioned,  the  assumed  deflection  configu- 
ration which  is  used  in  the  first  approximation  can  be  obtained 
by  multiplying  the  deflection  predicted  by  small-deflection 
theory  by  a magnification  factor  1/(1  - P/Pcr,)»  The  lowest 
critical  pressure  for  a circular  ring  under  uniform  external 
pressure  is  equ?l  to  3EI/R0^.  In  the  present  case  the  magnifica- 
tion factor  is  2.  After  the  bending  moments  and  the  forces  H 
and  V at  division  points  have  been  calculated  by  Marbec's  method 
[neglecting  the  terms  ( A0  )2  and  € ] the  angle-change  A 6 and 
the  deflections  ^ and  ^ of  small-deflection  theory  can  be 
obtained  by  means  of  Equations  (2. 2. 1.4)  and  (2.2. 1.9)  respectively. 
Multiplying  these  deflections  by  2,  a set  of  so-called  'magni- 
fied deflections'  are  obtained.  These  'magnified*  small-deflec- 
tions are  a very  good  approximation  to  the  true  or  large  deflec- 
tions. 

The  numerical  data  involved  in  this  method  of  analysis 
and  comparisons  of  the  results  obtained  by  large-deflection  and 
small-deflection  theories  are  tabulated  in  Table  (2. 2. 2.1). 

From  the  analysis  of  this  problem,  it  can  be  seen  that 
the  effect  of  the  extensibility  of  the  ring  axis  is  extremely 
small. 


(2) . Ring  Having  Shape  of  Third  Buckling  Mode 
A ring  having  an  irregularity  corresponding  to  the  shape 
of  the  third  buckling  mode  of  a circular  tube  is  given  by  the 
expression 
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f = Ro(l  - <f cos6p  ) , 

The  shape  Is  shown  In  Fig.  2.2.1. 5*  Two  cases  £ — 0.01  and 
£ **•  0.1  are  investigated.  Since  the  ring  is  symmetrical  with 
respect  to  both  X and  Y axes,  only  one  quadrant  need  be  investi- 
gated. Six  equal  angular  divisions  are  used,  as  in  previous 
examples. 

The  extension  of  the  ring  axis  is  assumed  to  be  negligi- 
ble in  this  example.  The  governing  equation,  then,  becomes 

MoJ'gj'+Vo 5 + /ws ^ t = 0 . 

The  curved  length  factor,  B — ds/d^>  , is  given  by 

B*R0(1  + 36  2 $ eos6  - 35  (S^cos 26?  )*. 

The  slope  of  the  tangent  to  the  ring  axis  at  any  point  is  ex- 
pressed by  the  formula 

tan$=(l  -£cos6<p+  6 £ sin6  f>  tan6lp)/[6  £ sin6<p  -(l-£cos6^)tan6p] . 

For  S =■  0.01,  the  numerical  data  are  tabulated  in  Table 
(2.2.2.2a) . The  final  results  for  bending  moment  and  defleotion 
for  both  cases  £=0.01  and  £=0.1  are  shown  in  Table  (2.2.2.2b). 

(3) . Circular  Ring  Having  Two  Symmetrical  Dents 

The  circular  ring  having  two  symmetrical  dents  examined 
in  Article  (2.1.2)  is  re-analyzed  using  large-defleotion  theory 
in  this  article.  Cases  15*,  £ =f  o,01  and  f^**15*, 

«£  = 0.1  under  external  and  internal  pressure  are  investigated. 

It  is  assumed  that  the  extension  of  the  ring  axis  is  negligible. 
The  slope  of  the  tangent  to  the  ring  axis  at  any  point  is 
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given  by: 

for  lP<lPo  tane=  Cf + tan  <*>)/(  -()-Un(p) 

and  (J)  > *Po  tan  © = — x/y 

Due  to  symmetry  only  six  equal  angular  divisions  in  one 
quadrant  of  the  ring  (with  ten  subdivisions  .in  the  dented  por- 
tion) need  be  investigated. 

If  the  pressure  acting  on  the  ring  is  internal,  the  govern- 
ing equations  oan  be  written  as  follows: 

H0  =»  0,  H6  =»  p(y*  ♦ 7)i  ) 

V0  « -p(x„  f ),  V6«0 

and 

For  rings  such  that  the  depth  of  the  indentation  is  one- 
hundredth  of  the  nominal  radius,  <T=s  0.01,  subjected  to  external 
and  internal  pressure,  numerical  data  are  given  in  Tables 
(2. 2. 2. 3)  to  (2. 2. 2. 5)  and  the  bending  moments  and  deflections 
for  all  cases  are  tabulated  in  Tables  (2.2.2.4b)  and  (2.2.2.5b). 
The  bending  moments  and  deflections  of  a ring  under  internal 
pressure  in  small-deflection  theory  are  not  given  in  the  table. 

For  small-deflection  theory  the  distinction  between  internal  and 
external  pressure  merely  creates  a change  in  sign. 

It  is  quite  interesting  to  note  that  the  influence  of 
deflection  of  the  ring  axis  on  the  final  bending  moment  and 
deflection  is  different  for  rings  having  various  kinds  of  imper- 
fections. For  an  elliptical  ring  having  an  ellipticity  equal  to 
0.01,  both  the  bending  moments  and  deflections  at  division  points 


to 


increase  100$  compared  with  those  given  by  small-deflection 
theory  for  p /per.  — 0.5.  On  the  other  hand,  in  the  case  of  a 
ring  whose  shape  corresponds  to  the  third  buckling  mode,  the 
maximum  bending  moment  is  only  a few  percent  greater  than  that 
obtained  by  Marbec's  method.  The  deflection  component,  either 
horizontal  or  vertical,  calculated  from  large-deflection  equations 
also  is  unaffected.  For  the  ring  having  two  symmetrical  dents, 
the  maximum  bending  moment  in  the  ring  increases  17$  for  £=-0.01 
and  10$  for  S =s  0.1.  However,  the  horizontal  deflection  increases 
about  75$  and  the  vertical  deflection  is  about  110$  greater  than 
that  predicated  by  small-deflection  theory. 

If  the  ring  is  acted  on  by  internal  pressure  instead  of 
external  pressure,  the  axial  tensile  force  developed  in  the  ring 
tends  to  reduce  the  deflection  produced  by  the  Internal  pressure. 
Therefore,  consideration  of  the  deflection  of  the  ring  axis  in 
the  large-deflection  theory  will  result  in  bending  moments  and 
deflections  less  than  those  predicted  by  small-deflection  theory. 
In  the  analysis  of  the  ring  having  two  symmetrical  dents,  for 
oases  $ =*  0.01  and  $ = 0.1,  the  horizontal  and  vertical  deflec- 
tions are  about  25$  to  30$  and  35$  to  40$  respectively  less  than 
those  in  the  small-deflection  case.  The  maximum  bending  moment 
drops  off  about  6$  for  S =*  0.01  and  12$  for  J s=  0.1. 

From  the  numerical  data  of  illustrative  problems  analy- 
zed by  the  large-deflection  equations,  besides  noting  that  the 
extension  of  the  ring  axis  is  extremely  small  and  can  be  neglected 
in  the  analysis,  it  is  also  worthwhile  to  point  out  that  the 
ohanges  of  lines  of  application  of  uniformly  distributed  external 


loads  (due  to  the  defleotion  of  the  ring  axis)  do  not  materially 

affeot  the  magnitude  of  the  bending  moment.  The  effect  of  axial 

foroes  in  the  ring  is  the  only  essential  one  which  does  play  a 

significant  part  in  the  final  determination  of  bending  moment 

and  defleotion.  That  effeot  is  great 

when  (a)  p /p  is  large,  say  greater  than  0.2. 
or « 

and  (b)  The  irregularity  resembles  the  oval  shape  of  the 
lowest  bucking  mode  of  a circular  ring. 
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2.3.  Simplified  Large-Deflection  Theory 
2.3.1.  Method  of  Analysis 

Due  to  the  fact  that  the  numerical  method  devised  for  the 
large-deflection  analysis  of  a circular  ring  having  initial 
irregularities  requires  a considerable  amount  of  numerical  work, 
it  has  seemed  neoessary  to  develop  a simplified  method  by  means 
of  which  a rapid  analysis  of  such  problems  can  be  made.  In  the 
case  of  a circular  ring  having  initial  out-of-roundness  under 
the  action  of  uniform  normal  pressure,  either  external  or  inter- 
nal, a simplification  is  effected  by  using  the  bending  moments 
predicted  by  the  small-defleot ion  theory  and  multiplying  them 
by  factors  proportional  to  the  relative  amplitudes  of  the  var- 
ious buokling  modes  present  in  the  original  imperfect  shape. 

The  maximum  bending  moments  calculated  by  this  rapid  method  of 
analysis  agree  quite  satisfactorily  with  the  results  obtained 
from  the  numerical  solution  of  the  large-deflection  equations. 

The  basic  concept  of  the  method  of  analysis  presented  in 
this  article  is  not  entirely  new.  From  the  elementary  theory 
of  bending  it  is  well  known  that  if  a beam  is  submitted  to  the 
action  of  lateral  loads  alone  the  deflection  of  the  beam  and  the 
stress  produced  are  proportional  to  the  magnitudes  of  the  loads. 
Small  changes  in  the  positions  of  the  loads  have  only  a small 
effect  on  deflections  and  stresses.  Small  deflections  do  not 
materially  affect  the  bending  moment  in  the  beam  and  in  the  cal- 
culation of  these  quantities  it  is  the  usual  practice  to  neglect 
entirely  the  deflections  of  the  beam  resulting  from  the  action 
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of  external  loads. 

In  case  axial  and  lateral  loads  act  simultaneously,  con- 
ditions are  entirely  different.  The  stresses  and  deflections 
are  no  longer  proportional  to  the  magnitude  of  the  longitudinal 
force.  A slight  eccentricity  in  the  application  of  the  axial 
load  or  a slight  deviation  of  the  axis  of  the  beam  from  a straight 
form  may  have  a substantial  effect  on  the  deflections  of  the 
beam  and  on  the  stresses  in  it.  The  effect  of  small  lateral 
deflections  of  the  beam  is  no  longer  negligible  and  must  be  con- 
sidered in  the  calculation  of  bending  moments.  When  the  magni- 
tude of  the  axial  compressive  force  approaches  a certain  limit- 
ing value,  usually  designated  the  ’critical’  load,  the  deflection 
becomes  very  sensitive  to  small  changes  in  the  position  of  the 
point  of  application  or  in  the  magnitude  of  the  axial  load. 

When  a straight  bar  is  under  the  simultaneous  action  of 
lateral  and  axial  loads,  the  maximum  bending  moment  or  deflection 
is  correctly  obtained  by  using  the  bending  moment  or  deflection 
produced  by  lateral  loads  alone  multiplied  by  a certain  factor. 
The  value  of  this  factor  depends  only  on  the  magnitude  of  the 
ratio  of  axial  load  to  its  critical  value.  It  approaches  unity 
as  the  compressive  force  becomes  smaller  and  smaller.  It 
Increases  indefinitely  when  the  compressive  force  approaches  the 
critical  value.  Without  introducing  serious  error,  the  factor 
can  be  replaced  by  a well-icnown  simple  form,  1/(1  - p/pCr.)» 

A ring  subjected  to  the  action  of  uniform  pressure  cor- 


responds to  the  case  of  a curved  beam  under  the  simultaneous 


action  of  lateral  and  axial  loads  the  latter  of  whioh  is  unknown 
in  magnitude  and  varies  along  the  ring  axis.  Although  the  ring 
problem  is  more  oomplioated  than  the  oase  of  a straight  bar,  the 
fundamental  idea  is  still  valid.  Due  to  the  presenoe  of  axial 
forces  in  the  ring,  small  deflections  of  the  ring  axis,  resulting 
from  the  action  of  external  pressure,  have  an  important  effect 
on  the  final  deflections  and  bending  moments.  We  can  say  that 
the  maximum  bending  moment  or  maximum  deflection  which  is  pro- 
duced by  the  action  of  uniform  pressure,  can  be  obtained  by 
multiplying  the  bending  moment  or  deflection  obtained  from  small- 
deflection  theory  by  certain  factors.  The  problem  is  how  to 
determine  the  magnitudes  of  these  factors  and  the  way  they  are 
to  be  used. 

Any  dosed  curve  defined  by  ^=R0[l  - f(<P)]  can  be 
expanded  in  a Fourier  series  (19)  of  buckling  mode  shapes,  as 
follows: 


where 


? = - z»  (A  cosn<f+  br>  s*nn<p)] 

n = i,2, — 

^ = 15  [.I 

an  - j(<p)Ctsnrif 

b„  = -yj* 


If  the  curve  is  symmetrical  with  respect  to  one  coordinate  axis, 
Equation  (2. 3.1.1)  becomes 
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In  the  ease  of  a curve  having  two  axes  of  symmetry,  the  odd  modes 
in  Equation  (2. 3.1.2)  disappear. 

This  equation  implies  that  a ring  of  any  shape  is  equiva- 
lent to  a sum  of  rings  having  the  shapes  of  buokling  modes.  The 
bending  moment  for  eaeh  mode  shape  can  readily  be  determined 
by  Marbec’s  method  and  the  sum  of  the  bending  moments  will  be 
the  total  bending  moment  developed  in  the  given  ring.  Consider- 
ing the  effect  of  deflection  of  the  ring  axis  resulting  from 
external  pressure,  the  maximum  bending  moment  developed  in -the 
ring  can  be  obtained  by  summing  up  all  the  bending  moments  of 
the  buokling  mode  shapes  present  in  the  original  imperfeot  form 
multiplied  by  their  proper  magnification  factors. 

The  critical  pressure  for  a circular  ring  under  uniform 
pressure,  for  each  buckling  mode  shape,  can  be  found  by  the 
formula : 

Pcr#=(n2  - 1)EI/R03  for  n ~2,  3,  4, 2. 3. 1.3 

The  magnification  factor  for  each  buckling  mode  shape  is  there- 
fore expressed  by 

F(n)  =.i/(i  - p/pcr#)  =*[l  - ^/(n2  - 1)]"'  2. 3.1.4 

where  ^ =*  pR0^/EI. 

The  magnification  factor  given  by  Equation  (2. 3. 1.4)  is  to  be 
used  in  the  simplified  method. 

For  a given  value  of  \ the  magnification  faotor 
approaches  unity  for  higher  values  of  n,  the  mode  order.  For  most 
practical  cases,  the  magnification  factors  of  a few  of  the  lower 
order  components  are  of  importance;  the  others  can  be  taken  as 
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unity.  The  problem,  then,  is  reduced  to  the  determination  of 
the  bending  moment  predicted  by  small-deflection  theory  and  a 
few  components  of  the  bending  moment  for  the  lower  mode  shapes. 

If  the  ring  has  a shape  which  corresponds  to  a pure  buck- 
ling mode,  the  bending  moment  in  the  ring  involves  only  the  one 
component  of  that  mode.  The  maximum  bending  moment  predicted  by 
this  rapid  method  of  analysis  is  obtained  by  using  the  bending 
moment  for  the  small-deflection  case  multiplied  by  one  magnifica- 
tion factor  corresponding  to  the  particular  buokling  mode  in 
question.  If  this  mode  is  a higher  one,  the  factor  is  approxi- 
mately equal  to  unity;  therefore,  the  maximum  bending  moment  pre- 
dicted by  this  method  will  differ  very  slightly  from  the  bending 
moment  predicted  by  small-deflection  theory.  This  was  observed 
to  be  the  oase  for  the  ring  treated  as  example  (2)  in  Article 

(2.2.1)  (see  Table  2. 2. 2. 2). 

If  the  initial  irregularity  of  a ring  is  very  small,  a 
further  simplification  can  be  made.  As  mentioned  in  Article 

(2.1.1) ,  the  bending  moment  at  any  point  on  a ring  is  equal  to 
a constant  force  (the  axial  force  in  a perfect  circular  ring) 
multiplied  by  the  amplitude  of  the  irregularity  at  that  point 
measured  with  respect  to  the  node  circle.  The  node  circle  is 
auoh  that  the  areas,  on  either  side  of  the  ring,  bounded  by  the 
ring  axis  and  the  node  oirole  balance  each  other.  Now  in  the 
Fourier  expansion,  the  quantity  R0(l  - a0/2)  has  a physloal  sig- 
nificance. The  t'erm  a0/2  which  is  (1/21T)  f*  f ( y )d«p  represents 
a correction  which  will  equalize  the  positive  and  negative  areas 
bounded  by  the  given  curve  t((f)  and  a reference  curve.  In  other 
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words,  the  radius  of  the  node  circle  is  R„(l  - ar/2). 

The  expression  an  C08fc^) » which  is  equal  to 
R„(l  - a0/2)  - £ , also  has  a physical  meaning.  The  coefficient 
of  the  expression,  an  , for  different  values  of  n will  give  the 
amplitude  of  the  imperfection  as  far  as  the 'corresponding  buck- 
ling mode  shape  is  concerned.  This  quantity,  for  small  imper- 
fections, is  proportional  to  the  magnitude  of  the  bending 
moment.  In  view  of  this  simplification,  even  the  small-deflec- 
tion analysis  is  unnecessary  and  the  only  work  required  is  that 
of  expanding  the  given  shape  of  the  ring  into  a Fourier  series. 

The  coefficient  of  the  zero’ th  mode  then  defines  the  radius  of 
the  node  circle  and  the  coefficient  of  the  n’th  mode  determines 
the  bending  moment  in  that  mode. 

The  technique  mentioned  in  this  article  is  also  applicable 
to  the  case  of  a ring  subjected  to  the  Internal  pressure.  The 
’magnification’  factor,  which  is  1/(1  + p/p„_  ) in  this  case, 
is  less  than  unity. 

2.3.2.  Illustrative  Problems 

In  order  to  illustrate  the  application  of  the  simplified 
method  of  analysis  developed  in  Article  (2.3.1)  to  a practical 
problem,  the  same  numerical  examples  given  in  Artiole  (2.2.2) 
are  reworked. 

(1).  Elliptical  Ring  (Fig.  2. 1.2.1) 

The  ellipse  (x/0.99Ro)2  ♦ (y/l.01Ro)2  = 1 differs  very 
slightly  from  the  lowest  buckling  mode  shape,  ^—R0(l  - 0.01oos2y>). 
The  bending  moment  in  the  elliptical  ring  is  the  sum  of  moments 
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for  each  of  the  mode  shapes.  The  lowest,  of  course,  predomin- 
ates because  it  has  both  the  largest  an  and  the  largest  magni- 
fication factor.  For  simplicity,  all  other  modes  can  be 
negleoted  without  introducing  appreciable  error.  The  bending 
moment  at  point  6 predicted  by  the  rapid  method  of  analysis 
developed  in  this  dissertation  is  simply  the  bending  moment  pre- 
dicted by  small-deflection  theory  multiplied  by  the  magnifica- 
tion factor  for  the  lowest  mode.  In  other  symbols, 

M6  * -.010050  R0xpR«,x  1/(1  - 1/2}  ss  -.020100  pRQ2 

2 

This  may  be  compared  with  a value  of  -.019829  pRQ  obtained  by 
the  ’exact'  large-deflection  analysis.  Since,  however,  the 
initial  imperfection  of  the  ring  is  very  small  the  bending 
moment  for  small-deflection  theory  need  not  even  be  calculated. 

It  may  be  considered  that  the  expression  p = R0(l  - 0.01cos2p) 
describes  the  ellipse  (x/0.99Ro)  +•  (y/1.01Ro)  =»  1.  Comparing 

two  expressions 

p - Ro(0- y)- 005  n?]  and  f = R*(  l-O.oi  cos  Z*?) 

it  is  evident  that  a0=a  0,  a^rsO.Ol,  and  all  other  4n  = 0 • 

It  may  be  seen  that  the  radius  of  the  node  circle  is  Re  and  the 
irregularity  at  point  6 measured  with  respect  to  the  node  circle 
is  equal  to  -0.01Ro  . The  axial  force  in  a perfect  circular 
ring  is  pR0  so  that  the  approximate  maximum  bending  moment  at 
point  6 can  be  calculated  as  follows: 

* - 


0.01  R0*pR.Xl/(l  - 1/2)=  - 0.0200  pR* 
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(2) .  Ring  Having  Shape  of  Third  Buckling  Mode  (Fig.  2.2.1. 5) 
Since  the  ring  has  a pure  buckling  mode  shape,  the  bend- 
ing moment  developed  in  this  ring  contains  only  a single  compo- 
nent. Using  Marbec's  method  for  the  small-deflection  analysis 
the  following  results  are  obtained: 

For  5 = 0.01,  r « 1.000075  Rn,  M - 0.010025  pR  2 

and  for  5=0.1,  rc  =r  1.006901  R0,  MQ  = 0.101925  pR02 

The  magnification  factor  for  the  third  buckling  mode  shape  is 
1/(1  - 1 • 5/ 35)  =*  70/67.  The  maximum  bending  moment  (at  point  0) 
predicted  by  the  simplified  method  is 

for  5=0.01,  M0  =:  0.010025  R0XpR0  X 70/67  =>  0.010474  pR  2 

and  for  5=0.1,  0.101925  R,XpR«  * 70/67  = 0.106488  pR#2 

Since  the  initial  irregularity  is  quite  small,  the  approx- 
imate method  can  be  applied.  We  do  not  even  have  to  know  the 
results  of  a small-deflection  analysis.  The  equation  of  the  ring 
is  at  once  in  the  form  of  a Fourier  series,  with  ae  equal  to 
zero.  The  irregularity  at  point  0 with  reference  to  the  node 

circle  equal  to  5r  . Then, 

0 

for  5=0.01,  M0=t  0.01  R0xPRo  x 70/67  = 0.010448  pR02 

and  for  5=0.1,  0.1  RcxpR0  X 70/67  = 0.104478  pR02 

In  this  oase  the  tedious  'exact'  analysis  (Table  2. 2. 2. 2)  yielded; 
For  5=0.01,  M0  = 0.010339  pRQ2 

and  for  5=0.1,  M0  =&  0.104280  pRQ2 

(3)  . Ring  Having  Two  Symmetrical  Dents 

The  shape  of  a ring  having  two  symmetrical  dents  is  given 
by  the  expression  (Fig.  2,1 .2  .2): 
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f = Ro[l-  S t05’(TTf/2^)] 

for 

-<£<<?<<£ 

e=R. 

for 

^2  ir-<Po 

p=  S'cV-ir/jr/afi] 

for 

7T-  ft  ^ 5 TT  f 

f = R. 

for 

7T+&  £ P < ZTT - to 

Due  to  the  symmetry  of  the  ring  with  respect  to  two  axes, 
the  Fourier  series  expanded  for  this  ring  is  given  by  Equation 
(2. 3.1. 5) 


0O 


f = Rolo-21)-  2 ‘*»c*n(p] 


2. 3. 1.5 
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The  Fourier  series  of  the  ring  is 
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2. 3.1.6 
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j 

(i).  and  0.01 

I 

By  substituting  $£*15°  e^d  J-^0.01  into  Equation 
(2. 3. 1.6) , we  have 


^ A * I *1  . 

P=?  R.[Cl-A000707)-  0.0127322  x jr-  CftsnyJ  2. 3.1.7 

Setting  n=2,  4,  - — in  Equation  (2. 3.1.7),  the  expression  for 
each  mode  shape  in  the  ring  is  readily  obtained.  For  instance, 

= Ro(0*999293  - 0.001396cos2  f ) 
fc4>  = Ro(0. 999293  - 0.001339cos4^  ) 

— Ro(0. 999293  - 0.001250cos6f  ) 

A single  problem  is  now  resolved  into  a set  of  problems, 
each  of  which  is  independent  of  the  others  and  can  be  solved 
separately.  The  bending  moments,  at  point  0,  for  the  first 
three  mode  shapes  analyzed  by  Marbec’s  method  are  as  follows: 

Mq(2)  = 0.001394  pR02 

M ^ = 0.001338  pR  2 
o o 

= 0.001248  pR  2 

/ ^ j 2 

The  summation  of  MQ<  must  be  equal  to  .009246  pR#  , the  bending 
moment  at  point  0 of  the  ring,  which  has  been  oomputed  in  Article 
(2.1.2)  by  small -deflection  theory.  Hence, 

Mo8>+  Mc°0)+  Mol,+ 0.005266  pR*  , 

The  magnification  factors  for  these  modes  are 


I 
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FU)=2,  F(4>  = 10/9,  F(4)=  70/67  and  Fiay=  F°o)  = =1  . 

Since,  for  higher  mode  shapes,  the  magnification  factor  approaches 
unity,  the  bending  moment  at  point  0 predicted  by  the  simplified 
method  will  be 

M*  =s  (.001394  X 2 4 .001338  X10/9  + .001248X7C/67  4 0.005266)pl^ 

= (.002788  4 .001486  4 .001304  4 .005266  ) pR* 

= .010844-  pR  2 

e 

If  only  the  first  component  of  bending  moment  were  considered, 
the  maximum  bending  moment  at  point  0 would  be  computed  as 
Mo  =•  (.001394  X 2 4 .007852)  pR02 

= 0.010640  pR02  . 

The  discrepancy  between  these  two  results  is  about  2%. 

Since  $ is  very  small,  a further  simplification  can  be 
made.  This  device,  as  in  previous  examples,  makes  it  unnecessary 
to  know  the  results  of  a small-deflection  analysis  in  advanoe* 

For  ^*-0,  Equation  (2. 3. 1,7)  reduces  to 

f=R0[l  - 0.000707  - (0.001395  f 0.001339  4 0.001250  4-*)] 
= R„(0. 999293  - 0.009293). 

It  has  been  pointed  out  that  0.999293  R*  can  be  tsken  as  the 

radius  of  the  node  cirole  and  the  coefficients  0.001395  R«  , 

oan  be  considered  as  the  amplitudes  of  the  irregularity  for  cor- 
responding mode  shapes,  at  point  0,  measured  with  respect  to  the 

node  cirole.  The  approximate  bending  moment  at  point  0 in  the 

2 

small-deflection  oase  is  given  by  the  term  0.009293  pR#  . The 
approximate  maximum  bending  moment  at  point  0 analyzed  by  the 
simplified  method  is 


M0* (0.001395  X 2 * 0.001339  X 10/9  + 0.001250x70/67  + ---)pR02 

— (0.002790  + 0.001488  + 0.001306  f 0.005309  + ) pRQ2 

**  0.010893  pR  2 
o 

(11).  (pQ  =*  15®  and  $ 0.1 

In  the  case  =*  15#»  S » 0.1,  the  mode  shapes  of  the 
ring  are 

R„( 0.992926  - 0.013956cos2^  ) 
fW=rRo(0. 992926  - 0. 013394oos4^  ) 


The  bending  moments  at  point  0 for  the  first  two  mode  shapes 
(analyzed  by  Marbec’s  method)  are 

~ 0.013905  pR02 
Mi4)  = 0.013344  pR#2 

and  the  total  bending  moment  at  point  0 which  was  computed  in 
Article  (2.1.2)  for  the  small -deflect! on  case  is  equal  to 
0.087926  pR02. 

Then,  the  maximum  bending  moment  at  point  0 is 
Mc  3.  (0.013905  X 2 + 0.013344  X 10/9  + 0.060677)  PR0 
— 0.103313  pRQ2  . 

When  the  ring  is  acted  on  by  internal  pressure,  a negative 
value  of  p must  be  used  in  computations.  The  magnification 
factor,  1/(1  + P/Por#)  is  loss  than  unity. 

For  y£e»15*  and  S ■=*  0.01,  the  maximum  bending  moment  at 
point  0 predicted  by  the  simplified  method  is 

M.  = -(0.001394  X 2/3  f 0.001338  X10/11  + 0.001248  X70/73 

0.005266)  pRo2 

= -.008608  pR  2 
0 * 


In  the  case  (p0  =*  15*  * $ =*  0.1,  the  maximum  bending  moment 

is  equal  to 

Mo  - -(0.013905  X 2/3  4 0.013344 x 10/11  4 0.060677)  pR  2 

o 

0.082080  pR  2, 
o 

The  comparison  of  results  for  the  illustrative  problems 
analyzed  by  the  different  methods  given  in  this  dissertation 
are  tabulated  in  Table  (2. 3.2.1). 


3.  CYLINDRICAL  SHELL  HAVING  SMALL  INITIAL  IRREGULARITIES 


3.1.  Baslo  Assumptions 

The  method  introduced  in  the  proceeding  chapters  for 
the  analysis  of  a circular  ring  having  small  initial  irregulari- 
ties and  subjected  to  a uniform  pressure  can  also  be  applied  te 
the  case  of  a thin  cylindrical  shell  provided  that  the  initial 
imperfection  is  uniform  along  the  length  of  the  shell.  It  is 
only  necessary  to  substitute  the  flexural  rigidity,  D =* 
Elr/12(l  "V1),  for  SI.  A natural  question  arises  as  to  the 
limits  of  the  validity  of  this  procedure  when  applied  to  local- 
ized dents  in  cylindrical  shells.  This  question  has  been 
answered  (at  least  partially)  by  the  solution  of  the  shell 
equations  for  a dent  varying  in  magnitude  along  the  longitudinal 
axis  of  the  cylinder. 

The  middle  surfaoe  of  the  shell  under  consideration  is 
assumed  to  be  defined  by  the  equation: 

^ = R0[l  - $ (cos  «ix)(cos  nf)]  3.1.1 

where  & , m,  and  n are  parameters  which  determine  the  geometri- 
cal shape  of  the  shell. 

The  shell  deviates  from  its  original  circular  shape  by 
an  amount  <T oos{nx)  and  has  a oross  section  similar  to  the  n'th 
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buckling  mode  shape.  Its  variation  along  the  longitudinal  axis 
may  be  specified  by  the  wave  length  in  that  direction,  2L  *=■ 
2tt/j n.  When  the  wave  length,  2L,  is  large  compared  with  the 
mean  radius,  RQ,  of  the  cylindrical  shell  (L/R0»  1)  the  situa- 
tion approaches  the  two-dimensional  case  of  a cylindrical  shell 
having  a cross  section  similar  to  the  n*th  mode  shape  all  along 
its  length.  If  the  wave  length  is  small  (L/R0<<  1)  the  dent  in 
the  cylindrical  shell  is  highly  localized  and  the  shell  has  a 
corrugated  surface.  The  results  to  be  developed  are  based  on  a 

theory  whioh  does  not  hold  for  very  highly  localized  dents  so 

not 

that  results  mustAbe  pressed  to  this  limit.  When  the  parameter 
n is  zero  the  shell  has  a circular  cross-section  having  differ- 
ent diameters  at  different  seotions  along  the  length.  This  last 
case  has  been  considered  by  E.  P.  Burmistrov  (16). 

In  all  the  subsequent  discussions  the  following  basic 
assumptions  are  made: 

(1) .  The  initial  irregularity  of  ti®  shell  is  so  small 
that  second  and  higher  order  terms  in  a can  be  neglected. 

(2) .  The  square  of  the  quantity  (m^R»)  is  negligible  com- 
pared with  unity,  i.e.  the  wave  length  of  the  dent  in  the  longi- 
tudinal direction  is  large  in  comparison  with  the  mean  radius  of 
the  cylindrical  shell. 

(3) .  The  square  of  the  quantity  (n<£  ) must  also  be  small 
compared  with  unity,  i.e.  the  shell  must  have  the  cross  section 
of  one  of  the  lower  mode  shapes. 

(4) .  The  shell  is  a thin  shell;  its  thickness  is  so  small 
as  oompared  with  its  initial  radii  of  curvature,  R*  and  R^  , 
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that  the  quantities  h/Ry  and  h/R  oan  be  disregarded. 

(5) .  The  method  of  analysis  is  based  on  a small-deflection 
theory  in  which  extension  of  the  middle  surface  of  the  shell 

is  neglected. 

(6) .  The  normal  foroes  N*  , and  shear  must  be 

small  in  comparison  with  their  critical  values. 
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3.2.  Fundamental  Quantities 

As  is  conventional  in  developing  the  governing  equations 
for  a cylindrical  shell,  we  consider  a small  element  formed  by 
two  pairs  of  adjacent  planes  which  are  normal  to  the  middle 
surface  of  the  shell  and  which  contain  its  principal  curvatures. 
The  location  of  this  element  is  defined  by  the  coordinate  x and 
the  angle  . The  coordinate  axes  X and  Y tangent  at  0 to  the 
lines  of  principal  curvature  and  the  axis  Z normal  to  the  middle 
surface  are  shown  in  Fig.  (3.2.1).  The  curved  length  of  the 
Infinitesimal  element  is 

(ds)2  « A2(dx)2  4 B2(df )2.  3.2.1 

If  the  middle  surface  of  the  shell  is  given  by  ),  the 

symbols  A and  B represent  the  expressions: 

A«[l  4 (d  f /ds)2]2  and  B =.  4 (df  /&<p  )2]  * 3.2.2 

The  principal  curvatures  of  the  surface  are 

(l/Rx)=-^(d2f  /4*2)/*3B 

and  (l/Rp)  = [>*4  2(d£/d  f )2-  f (d2f /df2)]  /AB  . 

V 3.2.3 

In  the  case  of  the  shell  given  by  Equation  (3.1.1) 

(neglecting  the  second  and  higher  order  quantities  mJ>R0  and 

n£),  Equations  (3.2.2)  and  (3.2.3)  become 

Asl,  B **  R (1  - <£oos$u$cosn^>  ) , 

(l/Rx)«  -m2£  Rocosfcu0cos({i^)  , 3.2.4 

and  (1/R^)=[l  - (n2  - 1 ) $ c os(pj$c 09(51  f )] /R0 . 

Only  first-order  terms  in  the  (small)  deformation  £ have  been 

retained  in  the  foregoing  expressions.  It  is  later  shown  that 


67 


the  stresses  due  to  a dent  whose  cross-sectional  profile  is  not 
a perfect  circle  depend  on  these  terms.  On  the  other  hand,  in 
the  case  investigated  by  Burmistrov  (16),  where  the  oross- 
section  of  the  shell  is  always  a perfect  circle  and  only  the 
radius  varies,  the  second-order  terms  must  be  retained.  In  this 
case  the  first-order  terms  in  <T  cancel  and  without  the  second- 
order  terms  the  solution  reduces  to  that  for  a perfect  circular 
cylinder. 

In  order  to  analyze  the  equilibrium  of  the  element,  the 
stresses  acting  on  the  plane  faces  are  resolved  in  the  directions 
of  the  coordinate  axes  (Fig.  3.2. la).  The  normal  stress  is  taken 
as  positive  when  it  produces  tension.  The  positive  directions 
of  the  components  of  shearing  stress  are  taken  as  the  positive 
directions  of  the  coordinate  axes  if  a tensile  stress  on  the 
same  side  would  have  tbe  positive  direction  of  the  corresponding 
axis.  If  the  tensile  stress  has  a direction  opposite  to  the 
positive  axis,  the  positive  direction  of  the  shearing  stress 
components  should  be  reversed.  The  bending  moments  are  taken 
as  positive  when  they  tend  to  decrease  the  initial  curvatures. 

The  positive  directions  of  twisting  moment  are  indicated  in 
Fig.  (3.2.1b). 

The  resultant  normal  and  shearing  forces  (per  unit  length 
of  the  normal  section)  in  the  longitudinal  and  circumferential 
directions  shewn  in  Fig.  (3.2.1a)  are: 

rh/l 

N*  - <rx(i-«/Rf)dz 

•/-h/l 

. rVi 

Nr  <tv 
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rW% 

Nx^  ^ J l/«tf)dl 

rti/i  . 

NW  = |L  tittU-Z/Wdl 

J-"/t 

/v  rh/i  3.2.5 

Qx  =*J  Zxi  (\-2/Rf)dl 

-"a 

tfi  U~Z/Rx)dz 

Due  to  the  fact  that  the  curvatures  of  the  shell  In  the  longi- 
tudinal and  circumferential  directions  are  different,  the  shear- 
ing forces  and  N^x  are  generally  not  exactly  equal  to 

each  other,  although  it  is  still  true  that  Vxf  = Zpx  • 

The  bending  and  twisting  moments  (per  unit  length  of  the 
normal  sections)  are  given  by  the  expressions: 

f* 

Mx-J  v*z(i- z/R<p) dz  ) 

Since  the  thickness  of  a thin  shell  is  small  in  compari- 
son with  the  radii  of  the  curvature,  the  terms  z/Rx  end  z/R^ 
in  Equations  (3.2.5)  end  (3.2.6)  disappear  in  the  evaluation  of 
the  resultant  noxmal  forces,  shears  and  bending  moments.  Then 
is  equal  to  . 

Using  Hooke's  law  and  neglecting  the  terms  z/Rx  end 
z/R/  , the  relations  between  dlreot  stresses  and  strains  at  any 


My^J  <ty2  (l-Z/R^dz 
Mft- J (\~~z/R*)di 


3.2.6 
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point  are  as  follows: 

0*  +V6f  ) 

0>  = E 0~  <f»  ¥\t€*)  3.2.7 

Md  - *vx  = E(  i+v)"Vx^  /z 

Similarly,  the  elastio  relations  between  stresses  at  any 
point,  at  a distance  z from  the  middle  surface,  and  ohanges  of 
curvature  are  given  by: 

<r*  =-Ez(i-vir,C**+v;tW 

<J>  = - EZ(|- + V%x)  3.2.8 

and  Cxf  — Cpx  ■*  - Ez  ( 1 + v)_l  %*p 

Substituting  Equations  (3.2.7)  and  (3.2.8)  into  (3.2.5) 
and  (3.2.6)  respectively,  the  normal  forces  and  moments  reduce  to: 


N*  » 


MX  =r 

M (p  sst 

and 

where 

K = 

K(  €x  i V€p  ) 

K(  6x  ) 

K(1  - V)  t*f  /2 

-d(  fa  * v fa  ) 

-D(  %f  i V fa  ) 

™ D(1  - V )fa<p 

(i  -v*)/Eh 


3.2.9 
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3. 3*  Relations  Between  Stresses,  Strains  and  Displacements 


In  Investigating  the  relations  between  strains,  changes 
of  curvature,  and  displacements,  the  deformation  of  an  element 
must  be  considered.  In  Love's  "Mathematical  Theory  of  Elasti- 
city" (7),  relationships  of  this  type  are  developed.  Recently, 
expressions  for  these  relations,  in  which  the  changes  of  curva- 
ture differ  slightly  from  those  derived  by  Love,  have  been  pub- 
lished by  W.  R.  Osgood  and  J.  A.  Joseph  (8).  However,  with  the 
assumption  that  the  middle  surface  of  the  shell  is  inextensible, 
the  expressions  for  changes  of  curvature  derived  by  Love  and  by 
Osgood  and  Joseph  agree.  In  the  present  investigation,  expres- 
sions for  changes  of  curvature  in  the  form  derived  by  Osgood  and 
Joseph  will  be  used.  The  relations  between  strains  and  displace- 
ments are  given  by  following  expressions: 


n 

AS  s<p 


U) 

Kx 


= +7fF 

^ L 1 u 05 

€(P  8 Jf  + A0  Vr  Rf 

y V db  _ _u_  , i 

A ax  A8  ax  A0  + 8 

y __  u a f \ \ . / j i _■>  v a a . i . i _ a a 

- A )*  l Kxhf>  ay  + a1  a<p 

4ii-(  -Ulw  -U  J4L 
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J-(jL 

26  W l Rt)  r2B'Rx 


+ _L  JLOL  _ JL  iA.  iliL-  _L  ii.  2*L 

T AB  ^x>f  a*5  >X  AB2  <?x 


3.3.1 

On  substituting  the  known  quantities  A,  B,  1/R x and 
1/R^  given  by  Equations  (3.2.4)  and  their  derivatives  with 
respect  to  the  coordinates  x and  <p  into  Equations  (3.3.1), 
these  expressions  reduce  to: 

€rx  =■  -4im2R,£i>Cos»»iXccsn<p 

53  -jjr  + £a»vwxc©snf  tSwmsininy  Cos  nf 

/if  » -|J-  + CoswiXCosntp  - ^MVsi«m«Cnsn<p 


Xx  =*  + $ m*RoU  si'**  m* cos 


nc^ 


f £ $(n  w»X  us  nV>  fr  5 coimx  si«  a«f  [n(  t\l-  ] 

%xr  = 

♦s—t i&f) 

4*s  Si* *nX Cosnip ^mCVj-llV  - ^5-J 4.  Gas mx  5i«n<p 


3.3.2 
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Substituting  Equations  (3.3.2)  into  (3.2.9) , we  have 
+ ^VmWSmmx  c*4 n</  i 

^r  = ,<[iip'^+vix+  S (-a^+vmNO“«*"w<:«'’P 

+ $ mU  SmrnX  eeancp  4 J 

N*-KftM**£** ] 

+$m(mvf^l-v)  sinmxeasn^  + ] 

+ ( Vw>R«1'-  !)■£-  sm  mx  cos  n <f  + ] 

R+  ^ 

u =*  rj  \L*  JL*.  4.  ±^L  iv.  _ 1=3L  iil 

M*p  UL  Rb  Mf  TRo  ax 

+ ^mn'Ti  (l-v)u  cosmx  sihrt(^  + ] 

3.3.3 

The  dotted  lines  in  Equations  (3*3*3)  represent  terms,  whioh 
eventually  vanish  from  the  equations  of  equilibrium  because  they 
involve  either  (i)  the  product  of  S and  v,  (li)  the  products 
of  S and  the  derivatives  of  v and  w with  respect  to  x and  <f> 
or  (iii)  the  products  of  S and  the  derivative  of  u with 
respeot  to  (f  . 
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3.4.  Equations  of  Equilibrium 

In  the  general  theory  of  shells,  the  equations  of  equili- 
brium, which  are  formed  by  equating  to  zero  the  resultant  and 
the  resultant  moment  of  all  foroes  applied  to  an  infinitely  small 
element  of  a thin  shell,  have  been  derived  by  A.  E.  H.  Love  (7). 
Love's  equations  of  equilibrium  differ  somewhat  from  those 
derived  by  W.  R.  Osgood  and  J.  A.  Joseph  (6).  If,  however,  the 
middle  surface  of  the  shell  is  assumed  to  be  lnextenslble  and 
two  non-linear  terms  (small  in  the  case  under  study)  are  neglected, 
the  two  sets  of  equations  are  identical.  Due  to  the  extremely 
complicated  mathematical  form  and  non-linear  nature  of  the  equa- 
tions of  equilibrium  in  the  most  general  shell  theory,  assumptions 
have  been  made  to  simplify  the  problem.  The  most  important  one 
is  that  the  ibrces  Nx  , Nf  and  are  very  small  in  comparison 

with  their  critical  values.  Under  these  circumstances  their 
effects  on  bending  oan  be  neglected.  With  this  assumption  all 
non-linear  terms  in  the  equations  ofv equilibrium  which  contain 
the  products  of  the  resultant  forces  or  resultant  moments  and 
either  the  small  displacements  u,  v,  w or  their  derivatives 
vanish.  The  equations  of  equilibrium  are: 


7k 


ir  * fa  $tl M*i’  ~ w X?  ' Jt  $ "r  * AB  $ M» + Q*  “ 0 

i &■ + -k  &' M» - i2#* -$r%r*+ - Tfi  £*-«*=< 0 

3.4.1 

By  substituting  Equations  (3.2.4) , (3.3.3)  together  with  their 
derivatives  into  (3.4.1)  and  in  this  way  eliminating  Qx  and 
, three  equations  in  the  dependant  variables  uf  v and  w oan 
be  obtained. 


+ sin mx  c*sn<f  + cosmycosn^p-f — 0 

K(¥  -0 + ~h  +°( 


+ hi  ii.JiJivtt' 

• /to'  *»i  ixi9 s * * 


2C  )tx  2R**  )xif< 

+ 5oDsmxa)sn^K[-  -Vi*Fn]w+ 

— $$my*x  c®4W<p£'^  K+  D fmn £}L2lE&. — jjm *0 

i ik+i  Jr£  A,-±±L  ?**-) 

b*+  Ri  >x>*  ^ ^ If  y /£  *«*  ?77?  r£ 

+K^‘^-t+i^+z' 

+£c«SWUfeaWp^D[  C~^vn  1 ? ^ C*J  - t"  <0 

^]tKp^w+WPB.u. 
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+ *0 

3.4.2 

The  three  independent  equations  given  by  Equations  (3.4.2)  are 
sufficient  to  determine  the  three  unknown  displacement  compo- 
nents. A solution  whioh  satisfies  Equations  (3.4.2)  will  satis- 
fy the  conditions  of  equilibrium  and  compatibility.  From 
Equations  (3.4.2),  it  oan  be  seen  that  the  terms  whioh  do  not 
contain  parameters  £ , m and  n will  constitute  a set  of  equa- 
tions for  a perfect  circular  cylindrical  shell.  Terms  contain- 
ing £ , m and  n represent  corrections  due  to  the  imperfections 
of  the  shell. 

All  terms  represented  by  dotted  lines  in  Equations  (3.4.2) 
contain  either  (i),  the  products  of  £ and  v or  (il)>,  the 
products  of  £ and  derivatives  of  v and  w with  respeot  to  x and 
(p,  or  (iii),  the  products  of  S and  derivatives  of  u with 
respect  to  <p  • 
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3.5.  Method  of  Analysis 

As  mentioned  in  the  previous  article,  the  governing 
equations  for  a oylindrioal  shell  having  initial  irregularities 
consist  of  the  equations  for  a perfeot  circular  oylindrioal 
shell  plus  the  correction  terms  due  to  the  imperfections  of  the 
shell.  These  terms  contain  the  parameters  £ , m and  n.  Now  it 
is  required  to  find  displacement  components  whioh  satisfy  the 
governing  equations  and  at  the  same  time  satisfy  the  boundary 
conditions  of  the  given  problem.  The  nature  of  the  problem 
suggests  trying  a solution  of  the  form 
u = U0  (x,  <(>  ) +$  U^x,^  ) 

v=  V0(x,ip  ) + $ V^x,^)  3.5.1 

w = W0  (x, (f  ) t $ Wx(x,p  }. 

Here  U0(x,p),  V*(x,^)  and  Vf0(x,<f)  are  the  solution  for  a 
perfect  circular  cylindrical  shell  and  U^tx,^),  V^tx,^)  and 
V1-*’  ) represent  the  corrections  due  to  initial  imperfections 
of  the  shell. 

Substituting  Equations  (3.5.1)  into  (3.4.2)  two  sets  of 
terms  are  encountered,  one  proportional  to  £ and  one  independent 
of  £ . Each  must  vanish.  On  setting  them  equal  to  zero  two  sets 
of  linear  partial  differential  equations  are  obtained.  The 
first  group  is  given  by  Equations  (3.5.2).  They  contain  as  depen- 
dent variables  only  U©(x,^),  V0(x,qp)  and  W0(x,^}  which  depend 
on  the  given  external  loading  condition  specified  by  X* , Y* , and 
Z*. 

These  equations  refer  to  the  perfect  circular  cylindrical 
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shell* 


>x*  2fc>  + ^Rj  *f*  R*  /T  + K u 

c/Jb^  JlYk  + I±H  ilHs  + JL  £2k  _ l JSiHl£5+A  +-L^l 
r w 2x4  * 2R*  axi^>  RJ  ty*  rJ  >*Vy>  T f^4  w 


. l-v  *XV.  »-y  a*y«  , Ro7* y*  n 

T 2 2Xl  2R.  ?*Pf  D 1 “ U 


cfi  2-.ffitl^I.\_rl,‘i!]!4t2ii  , i /ff.  . ) j’W. 
h'RJ  *f  (C  R.  J l>  a*  + z ,x»*‘  + R*  >?♦  + RJ-  >f> 

I V)  j.  / 1 y)  J ffft  _ l-V  1 , R»*  7 f = 0 

+v^ui  vjft*  5^iJ+  D £ u 


3.5.2 


The  second  set  of  equations  contain  Uj(x,p),  V((x,^)  and 
W((x,p  ) as  dependent  variables.  They  depend  only  on  the  solu- 
tion of  the  first  set  which  give  U0(x,f),  V„(x,^)  and  W«(x,^). 
This  second  set  give  the  corrections  due  to  the  initial  irregu- 
larity of  the  shell. 

On  introducing  the  assumption  that  the  extension  of  the 
middle  surface  of  the  shell  is  negligibly  small,  (i.e.  6X=0, 

e?  * 0 and  2*^  « 0 in  Equation  3.3.1),  Equations  (3*5*2) 
reduce  to  those  given  by  S.  Timoshenko  (20). 

If,  now,  a system  of  loads  aoting  on  the  shell  is  symme- 
trical with  respect  to  the  plane  x=0  of  the  shell  and  symmetri- 
cal with  respeot  to  the  axis  0 in  any  plane  of  the  cross 

section,  the  radial  displacement  at  points  4 <f  and  •if  must 
be  equal  to  each  other  and  in  the  same  direction,  and  the  cir- 
cumferential displacement  at  all  points  on  the  symmetry  axis 
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must  vanish.  To  satisfy  these  boundary  conditions,  Equations 
(3.5.1)  can  be  taken  in  the  form: 

u ■=■  U0(x,<f  ) + <£  U,  (x)cosnfp 

v=rV0(x,<p)  + £ V|  (x)sinn^  3.5.3 

and  w=.W0(x,^)  + S W,  (x)cosn^p 

The  longitudinal  displacement  u is  chosen  as  a function  of  cosn^ 
so  that  the  variable  (f  in  the  second  group  of  linear  partial 
differential  equations  oan  be  separated.  Then,  the  second  group 
of  equations  reduces  to  a set  of  linear  ordinary  differential 
equations,  in  the  independent  variable  x. 

-VRoW,'  f Ro  V/  + Ro*  U;  - 

= sinmv-i; mx  , 

-nR0X/+n(F+nl-i)W1  + ^(F  + i)RezV,*  - n*F V,  - \ R.U/ 

= jF[ii(>»l-l)4V  nmlR0*  ]-2n(vi tosmx  - mv 

+ [FnmR0  + ^Cl-vJw3R0VTi(vwlf^-iJmR#JU6sir»mv 
- 1C Wi *"+  (2 (F+n4- - n(\-v) R z V; % n F Vj 
+ [VF-n*(l-V)]RolX 

= [F[-2(-»it-i)  -2vwlR*] +2wx(**i -i;  + V-/J  m'RfJ  Wv  c*  m* 

t[F(«lR^+vn*-v}R#+(zmtR*^v+rf-i;nVmlRo* 

- £ F*> R*  + (w>R^ - v+  nl-  v r»\)»n*Ro  t-  r)*(t; w^R^-OroR  J Uft  *’"»  rnX 


3.5.4 
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When  the  external  loads  have  been  speoified,  Equations 
(3.5.2)  oan  be  solved  and  the  solution  U0(x,^},  V,[x,<p)  and 
W0(x,^)  used  for  solving  Equations  (3.5.4).  However,  the 
application  of  the  general  theory  of  bending  of  circular  cylin- 
drical shells  in  even  the  simplest  cases  may  result  in  very  com- 
plicated calculations.  To  make  the  theory  applicable  to  the 
solution  of  practical  problems  some  further  simplifications  in 
this  theory  are  necessary. 

The  shell  concerned  is  under  the  action  of  uniform  exter- 
nal pressure.  In  equations  (3.5.2),  X*  and  Y*  are  equal  to  zero 
and  Z’  is  equal  to  a constant,  p.  The  solution  for  the  case  of 
a circular  cylindrical  shell  subjected  to  uniform  external  pres- 
sure is  well  known. 

When  the  shell  is  free  to  expand  in  a longitudinal 
direction  at  both  ends  we  have  as  a solution  of  Equations  (3.5.2) 
• 

U„(x,<?)  * (tfpR«,/Eh)x  , V,  (x,  f ) * 0 , VMx,^)  = pR0/Eh 

When  the  shell  is  restrained  at  its  ends  the  approximate  solu- 
tion of  Equations  (3*5*2)  is 

Uo  a 0 , V,  a 0 , W«  a w0  ■ (1  -V*)pR*/Eh  . 3. 5. 5 

The  above  solution  <is , of  course,  only  correct  at  points  dis- 
tant from  the  supports.  In  the  vicinity  of  the  supports  local- 
ized longitudinal  bending  occurs.  Since  interest  centers  in 
conditions  at  points  remoted  from  the  supports,  it  is  satisfac- 
tory to  use  Equations  (3.5.5)  in  Equations  (3.5.4)  for  the  deter- 
mination of  U|,Vj  , and  Wf  . For  a long  shell  the  longitudinal 
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beading  due  to  end  restraint  dies  out  rapidly  (20)* 

Assuming  now  that  U((x),  V(  (x)  and  W,(x)  have  the  form: 
Uj(x)  — rmn  sin£u$ 

Vi  (x)  = $mn  cosfnx)  3.5.6 

and  W,  (x)  = otmn  oos^ix) 

When  x =s  0,  U,  (x)  satisfies  the  boundary  condition  that,  due 
to  the  symmetry  of  the  shell  with  respect  to  the  plane  x*  0, 
u at  x =*  0 must  be  equal  to  zero.  Substituting  Equations 
(3.5.6)  into  (3.5.4)  and  using  Uc  *■  0,  Vo  **•  0 and  W0  =*  w»  , 
the  three  simultaneous  ordinary  differential  equations  reduce  to 
three  simultaneous  algebraic  equations. 

iJrn&cCmn  fmn  = [mJ  l^+mR.CnV-^Ju)®  / 

Tdfcn'ftf+CF  + *>x-  0> Un  - [mxRoX  (F  +0^  + ^F]  f 

,THaR-[F0  + v)_(J.v)]ymn 

=■  |p[n(nx-'l)-»-viimlR0x]  - 2n  (ft*-  l)Ju)c  > 

- [m4/?e*  + mxR0*(2n’'-v3f(F  + n*f-nx3]^mn  +n[mlg^(i-v)+*F]  pr..n 

F hi R©£v F-nx0"^)J  y»nn 

=^F[- 2(tix-i) - 2 v mxRol]  + 2nx(y\*  \)  + m'vi (2vnx-v-\ )j u)0 

3.5.7 

In  the  derivation  of  the  relations  between  strains  and 
displacements  and  of  the  equations  of  equilibrium,  the  extension 
of  the  middle  surface  of  the  shell  was  included.  It  can  be  seen 
that  in  Equations  (3.5.7)  the  terms  -2n(nz  - l)w,  and 
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[ 2n,^  (n*  - 1)  4-  m1R  (2\?n  - \J  - 1)]  w0  in  the  right  hand 

sides  of  the  second  and  third  equations  represent  the  effect 
of  the  extension  of  middle  surface  of  the  shell.  In  the  case 
of  a shell  the  value  of  7 =»  12(R©/h)  is  very  large#  For  sim- 
plicity in  calculation,  the  terms  corresponding  to  the  effect 
of  the  extensibility  of  the  middle  surface  of  the  shell  can  be 
neglected. 

By  disregarding  these  terms  and  solving  Equations  (3.5.7), 
the  coefficients  o<mn  , jjwn  and  ifmfl  in  the  displace- 
ment components  are  obtained. 

*Un  = (G2/C1) (PR04/D) 

(*m«  = [nC3/(l  -1>  )C1‘]  (PRoVd)  3.5.8 

Ymn  --[mR^/U  jCjCpR^/D) 

where 

= F hx(*4v)4<> vVj*>uR* 

+ n%( ’4 3 )(w't- t) nfRo  4 n4 ( *1'- I )*f  p- ( maR*  4 ( -o) m * 
-vf(V>v-i)rn4R^  - 

C2~  4 3 t2^-2vV-Z+vJ»n4f?^ -(3-v)]***.1 

4 i ) 4-  I -V ) 4 v]mfcRj  4 l4iv  - V *■)  - 1 tvj  m4  R 4 

- nf  [2  *4(  J -v)  - y>Xx-  V ) + )]  j 

Cj*  (|-v)(|42i/)w4R04  4 (|-v)[2r?x(l4V)-j]rn*^4>ixCvil:-0(»“V) 
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+ “ - [ W rfii-v)  - ( i+  vV]  m*R06 

- [4n4(i-va)  -2oa(4-3vXl+»)t(i+v,}(5-5v  +2.v  *)]  m4R4 

- [2h6C|-v)(2+v)- n4(»i-2V  “7i)1)+w1'(ij-7»-2v^)-(5-5i/  +2Viy)Jr>i*'R* 

- Y)\l-vX*ifc-4*i*+5*»*-2)J 

+ -p  ^-O-vjwVf-  0 -*)[vnx~(  l«)]»»V>(l-»)[h4-n%(l-v^-v]  m4R4 
+ (_|-v)[vht-n4(l  ♦*)+«*■]  mlR*J 
C4.—  0 *vXl-  2Vl)m4R4  - (l  - vX  l-2v)m>&*+  nXT»v- 1 j 0 -*) 

+ -p  ^ C- v)m8R68+ (»- v)[anx(2fi0  - Cl  +v)]  w<,R,fc 
+ [4m4(I-v1)«h1'(4-|i>  -2v%)  + 0-v)(f+v-2i)x)J»»i4R4 
+ [2yo6(l-v)-n4(5’-tv-3v,9  +*x(5^9v)-(l-vXl'2v^m*R* 

- >lX[>l4(l-v)-2v,4t  >ix(l+3v;-2v]J 

- (1  - v)[n4-  n*0  - v*J  - v ]m4R?-(i-v)nx[v  n4-  rj4(  1 tv)  ♦ 1 J J 

When  the  deflected  surface  of  the  shell  under  the  action 


of  external  loads  has  been  determined,  the  internal  forces  and 
bending  moments  at  any  point  of  the  shell  can  readily  be  found 
by  means  of  Equations  (3.3.3)* 
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The  bending  moments,  Mx  and  , in  the  longitudinal  and 
circumferential  directions  respectively,  are  usually  the  deter- 
mining factors  in  the  failure  of  the  thin  shells.  These  bend- 
ing moments  at  any  point  on  the  shell,  have  the  following  ex- 
pressions respectively: 

Mx*  Casrup 

-251v(2n1-0^mn  Cos*mx  Cos xw<p 

+ £ [vn(vix-l)^»n  + ^hlolmr,  J CoS*W  Si'n’rupJ 
Mif  35  - '^^-^[(v'x~l)  + VtnxK0,'Jo(mn  c«amx  e«sn<p 
-ZSXln1^ 1)  dmn  Cos  vm*  Cas^np 
- 5 * [ ( I - v mvR  x) m R* }fmn  + m”  K.'oimn]  SmW  CoS* r»f 
+r[*K-  of  mn  fTlVwh]  CoSVroX  S«n Vn <pj 
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3.6.  Illixatratlve  Problems 


In  the  first  illustrative  problem,  the  surface  of  the 
cylindrical  shell  is  defined  by 

^ - S cos  (if  x/L)cos2<p]  . 3.6.1 

The  oross  section  of  the  shell  at  x«0  is  a first  buckling  mode 
shape  which  is  very  close  to  an  ellipse.  At  x*  L/2  the  imper- 
fection dies  out  and  the  cross-section  is  a circle.  At  a point 
one  half  wave-length  along  the  longitudinal  axis,  x=*L,  the  cross- 
section  of  the  shell  returns  to  its  elliptical  shape  with  a phase 
difference  of  90  degrees.  As  the  coordinate  x increases, the 
cross-section  approaches  a circle  again  and  when  x is  equal  to 
one  full  wave-length  the  cross-section  returns  to  its  elliptical 
shape  as  at  x-0. 

Neglecting  the  extension  of  the  middle  surface  of  the 
shell,  the  displacement  components  are  given  by  the  following 
expressions : 

u«  S s in  ( it  x/L)  cos2<p 

v=*  oos(  flr  x/L)sin2f  3.6.2 

ws=  i'daiz.  oos(  7T  x/L)cos2p 
and  the  values  of  Cs  in  Equations  (3.5.8)  are: 


C,  = m4 R4  F t njaR^  + l6m*K‘  +77^^  -H56maR*  + 144 
4 j-(rr)3R? t 4m*Re6  - 12  m4p  + - 48  m1  R*  ; 
!Om4R.4+2owi^+48  + y 
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Ci  +5**1?*  -M2  + ^n6Rftm*(£-57rn+Kf-i2JinxKZ-7z) 

+ (-  m8R*8+  m*  R?  1 1 2 m4  R? - \ 2 mXl) 

C4-  *n4R4-mxR.*  + 12  + -p-(*8R<?+  l5m6R0ft  +33  m4Rft4-6lwX1-)i^) 
+ -pO«8&8-»"6Ro 
where  m ■=*  tt/L. 

For  the  case  of  a thin  shell,  the  magnitude  of  F is  very  large 
and  the  terms  in  l/F2  can  be  neglected  in  the  calculation  of 
these  coefficients. 

The  bending  moments,  at  any  point,  in  the  longitudinal 
and  circumferential  directions  of  the  shell  are  given  by; 

Mx  =*  - -jg  (-5  •’^R.^cima  C<smx Coe2<p  Xm  sir\\nxcesxlp) 

Ccsnrrx Cos  " |4SV»(ni  coiXmx  Cosa2^> 

-5  (m  RoXmi  +ftiaRo  oiwi)  sm  vmx  cos*2y> 

+rc&p  ml  +4oCfnaJ  CoS^mX  Sin  1 2<f  J 

3.6.3 

The  maximum  bending  moments  in  the  longitudinal  and  circumferen- 
tial directions  oocur  at  the  point  x = 0 and  <p  =»  0.  They  are 
Max.  M*  <5*  m2fl<.miD 

3.6.4 

Max.  My,  =.  <T(3  + 14  J ) eU*D/Ro2 
in  which  only  the  coefficient  o(mz  appears. 
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The  maximum  bending  moments  Ux  and  My  at  xaO  and 
are  plotted  against  the  wave-length  of  the  dent  in  the  lon- 
gitudinal direction  for  R0/h  ranging  from  20  to  100  and  £ equal 
to  .01  in  Figs.  (3.6.1)  and  (3.6.2). 

For  values  of  Ro/h»60,  L/Ro=*10  and  S **  0.01,  the 
variation  of  bending  moments  and  along  the  line  <f  a*  0 
in  the  longitudinal  direction  of  the  shell  are  given  in  Figs. 
(3.6.3)  and  (3.6.4)  respectively. 

When  m is  equal  to  zero,  the  problem  reduces  to  a two- 
dimensional  one.  The  displacements  in  the  circumferential  and 
radial  directions  of  an  elliptical  ring  are  given  by. 

v s=<TpR0^sin2^/6D 
and  w « £pR0^cos2j?  /3D 

The  maximum  bending  moment  at  point  of  ^=3  0 is  equal  to 

Max.  My  — • (3  + 14  & ) <5  PR0V3  . 

In  the  case  0.01,  My  = 0.010467  pRQ^  which  differs  from  the 
result  obtained  by  Marbec's  method  by  about  5%. 

In  the  other  extreme  case  m is  infinitely  large,  i.e.,  the 
wave-length  of  the  dent  along  the  longitudinal  axis  of  the  shell 
approaches  zero.  In  this  case,  the  bending  moments  M * and  1 
given  by  Equations  (3.6.3)  tend  to  become  constant.  However, 
when  the  dent  is  highly  localized,  the  results  of  the  solution 
given  in  the  proceeding  article  are  not  reliable  sinoe  the  method 
of  analysis  is  based  on  the  assumption  that  the  wave-length  of 
the  dent  along  the  longitudinal  axis  of  the  shell  is  large  com- 
pared with  the  mean  radius  of  the  shell. 
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Two  other  oases,  n=s2,  \>=l/3  and  n«4>V=il/3  are  also 
Illustrated.  The  maximum  bending  moments  and  at  x = 0 

and  ^ « 0 are  similarly  plotted  against  the  wave-length  of  the 
imperfection  for  R«/h  ranging  from  20  to  100  and  & = 0.01.  By 
specifying  R©/  h«60,  L/R0  =*10  and  J"-=0.01  as  in  the  case  of 
n =»  2,  V =*  0,  the  bending  moments  and  along  the  line 

^=.0  can  also  be  plotted.  These  curves  are  shown  in  Figs. 
(3.6.5)  to  (3.6.12). 

The  foregoing  curves  show  that,  when  the  wave-length  of 
the  dent  lies  within  a certain  range,  which  depends  on  the  value 
of  R„/h,  bending  moments  in  both  the  circumferential  and  the 
longitudinal  directions  are  negligibly  small  and  membrane  stres- 
ses predominate.  As  the  wave-length  exceeds  that  range,  the 
circumferential  bending  moment  increases  appreciably.  If  the 
wave-length  increases  still  further,  the  circumferential  bending 
moment  Increases  until  the  solution  is  identical  to  the  two- 
dimensional  case. 

On  the  other  hand,  when  the  wave-length  decreases,  bend- 
ing moments  increase  over  the  intermediate  wave-length  at  which 
membrane  stresses  predominate.  The  dent  becomes  more  localized 
in  character.  The  method  of  analysis,  however,  is  not  accurate 
in  the  vicinity  of  a highly  localized  dent. 

Other  things  being  equal,  smaller  values  of  R©/h  produce 
greater  bending  moments  in  shells  having  the  same  wave-length  for 
the  dent.  In  the  case  n«*4,  bending  moments  produced  are  larger 
than  those  for  n»*2  if  the  wave-lengths  of  the  dents  and  the 
values  of  R©/h  for  these  two  cases  are  the  same. 
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4.  CONCLUSIONS 


4.1.  Summary 

Investigations  of  circular  rings  with  initial  out-of- 
roundness  and  local  dents  are  presented  using  both  small- 
deflection  (approximate)  and  large-deflection  ('exact')  theo- 
ries. 

In  the  small-deflection  case,  a method  due  to  Marbec  is 
recommended  as  the  most  simple  and  convenient  one  for  the  pre- 
diction of  bending  moments  and  stresses  due  to  dents  and  initial 
out-of-roundness  in  rings  under  hydrostatic  pressure.  When  the 
initial  imperfection  in  a circular  ring  is  small,  the  bending 
moment  at  any  point  on  the  ring  can  be  obtained  by  multiplying 
the  amount  of  irregularity  at  that  point,  measured  with  refer- 
ence to  a 'node  circle'  defined  in  the  text,  by  the  constant 
force  (pR0)  in  the  ring. 

A numerical  approach  which  is  a combination  of  two  well- 
known  procedures,  is  developed  for  the  'exact'  or  large-deflec- 
tion analysis  of  circular  rings  with  initial  imperfections, 
taking  into  consideration  the  deflections  and  the  extension  of 
the  ring  axis.  A method  of  successive  approximations  is  used  in 
conjunction  with  the  numerical  procedure  of  integration  due  to 
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N.  M.  Newmark  (16).  Although  it  requires  a considerable  amount 
of  numerical  work,  this  procedure  has  the  merit  of  simplicity, 
wide  adaptibility  and  a straightforward  arrangement  of  com- 
putions.  The  results  of  illustrative  examples  reveal  that  for 
that  practical  problems  the  effect  of  the  extensibility  of  the 
ring  axis  is  negligible;  also,  changes  in  direction  of  the  exter- 
nal loads  due  to  the  deflection  of  the  ring  axis  do  not  mater- 
ially affect  the  final  stresses.  In  the  case  of  a thin  ring  with 
small  imperfections  they  can  be  disregarded  without  Introducing 
appreciable  error.  However,  due  to  the  smell  deviation  of  the 
ring  axis  from  its  original  shape,  the  presence  of  axial  forces 
contributes  substantially  to  the  stresses  produced  in  the  ring. 

A rapid  method  of  analysis  teking  into  account  the  effect 
of  axial  forces  is  developed.  By  means  of  this  procedure  the 
bending  moments  or  stresses  produced  in  rings  with  small  initial 
Irregularities,  under  hydrostatic  pressure,  can  easily  be  pre- 
dicted. This  simplification  is  obtained  by  multiplying  the  bend- 
ing moments  or  stresses  determined  by  small-deflection  theory, 
by  factors  proportional  to  the  relative  amplitudes  of  the  var- 
ious buckling  mode  shapes  present  in  the  original  dent.  The 
shape  of  a ring  with  Initial  out-of-ro undness  can  be  expressed 
as  a Fourier  series.  If  the  imperfection  of  the  ring  is  small, 
the  zero’ th  term  in  the  Fourier  series  gives  the  difference  of 
the  radii  of  the  circular  ring  and  the  node  circle  of  Barbee’s 
method.  The  coefficient  of  the  n'th  term  in  the  series  repre- 
sents the  amount  of  irregularity  of  the  n'th  buckling  mode  shape 
measured  with  reference  to  the  node  circle. 
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The  range  of  validity  of  the  simplified  procedure  has 
been  verified  by  the  foregoing  numerical  method,  in  case  the 
maximum  amount  of  initial  irregularity  of  a circular  ring 
ranges  from  1/100  to  1/10  of  its  original  radius,  results  ob- 
tained by  these  two  methods  check  quite  well.  If  the  initial 
irregularity  becomes  large,  appreciable  error  will  be  involved 
in  the  result  obtained  by  the  simplified  method.  When  the  max- 
imum amount  of  initial  imperfection  reaches  two  tenths  of  the 
radius,  the  maximum  bending  moment  predicted  by  the  rapid  method 
of  analysis  is  thirty  per  cent  greater  than  that  calculated  by 
the  numerical  method.  However,  for  most  practical  problems  the 
maximum  amount  of  initial  imperfection  will  hardly  be  more  than 
one  tenth  of  the  radius.  Furthermore,  if  the  imperfection  is 
very  large,  residual  stresses,  which  are  assumed  to  be  negligibly 
small  in  the  present  investigation,  might  be  of  prime  importance. 

A small  Initial  irregularity  greatly  reduces  the  strength 
of  a circular  ring.  In  small-deflection  theory,  when  the  maximum 
amount  of  imperfection  is  equal  to  1/100  of  the  radius,  the 
maximum  fiber  stress  due  to  bending  moment  is  from  four  to  seven 
times  the  hoop  stress,  depending  on  the  ratio  of  the  nominal 
radius  to  the  thickness.  The  foregoing  figures  are  for  values 
of  R„/h  ranging  from  50  to  100.  Taking  into  acoount  the 
deflection  of  the  ring  axis,  the  effect  of  such  an  irregularity 
on  the  strength  of  a circular  ring  under  uniform  external  pres- 
sure becomes  even  more  pronounced  in  most  practical  cases. 

The  bending  moments  and  henoe  the  stresses  developed  in  an 
elliptical  ring  (using  the  simplified  method)  are  obtained  by 
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multiplying  those  predicted  by  small-deflection  theory  by  a 
factor  1/(1  - P/PCr.)*  The  banding  moments  developed  in  a ring 
having  local  dents  depends  not  only  on  the  magnitude  of  the 
external  pressure  but  also  on  the  geometrical  shape  of  the 
dented  portion.  For  a particular  dented  shape  illustrated  in 
Fig,  (2. 1.2. 2),  the  maximum  bending  moment  was  found  to  be 
about  18%  greater  than  that  predicted  by  small-deflection  theory. 
It  may  also  be  concluded  that  if  the  shape  of  a ring  is  similar 
to  that  of  the  lower  buckling  modes  the  effect  of  the  deflection 
of  the  ring  axis  on  the  stresses  will  be  of  primary  significance. 
For  example,  the  worst  case  examined  is  that  of  a ring  having  an 
initial  ellipticity  under  hydrostatic  pressure. 

Since  the  foregoing  solution  is  two-dimensional  in 
nature,  a question  arises  as  to  the  limits  of  its  validity  when 
applied  to  localized  dents  in  cylindrical  shells.  This  question 
has  been  answered  (at  least  partially)  by  the  solution  of  the 
shell  equations  for  a dent  varying  in  magnitude  along  the  longi- 
tudinal axis  of  the  shell.  Due  to  the  complexity  of  the  general 
shell  equations,  some  simplifications  have  been  made  in  order  to 
make  the  problem  solvable.  The  results  of  analysis  indicate 
that  when  the  wave-length  of  the  dent  is  large  compared  with  the 
initial  radius  of  the  cylindrical  shell,  the  solution  coincides 
with  the  two-dimensional  case.  As  the  wave-length  decreases, 
bending  moments  decrease  and  membrane  stresses  predominate.  As 
the  wave-length  decreases  still  further  and  the  dent  becomes  more 
localized  in  character,  bending  moments  again  Increase.  The 
method  of  analysis,  presented  in  this  dissertation,  however,  is 


not  accurate  in  the  vioinity  of  a highly  localized  dent. 
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APPENDIX  I 

Nomenclature 

The  following  nomenclature,  used  in  this  report, 
conforms  as  nearly  as  possible  to  the  conventional  symbols  of 
applied  mechanics.  The  quantity  in  parentheses  represents  the 
dimension  units  corresponding  to  the  symbol. 

E =•  Young's  modulus  of  elasticity  ( MT’2  If*  ). 

A =*  Area  of  cross  section  ( L2  } . 
h = Thickness  of  either  a ring  or  a shell  ( L ). 

Rq=  Radius  of  a circular  ring  or  a circular  cylindrical  shell  ( L). 
x,  y *ar  Coordinates  of  any  point  on  the  ring  axis  with  respect  to 
the  principal  axes  of  the  ring  ( L,  L ). 
x' , y'=.  Coordinates  of  any  point  on  the  ring  axis  with  respect 
to  any  arbitrary  axes  X'  and  Y*  (L,  L ). 
a , b =■  Coordinates  of  the  center  of  the  node  circle  with  respect 
to  the  principal  axes  { L,  L ). 

f>}(f  ■=  Coordinates  of  any  point  in  polar  coordinates  ( L,  Rad.  ). 
rc  Radius  of  the  node  circle  ( L )• 

pc  =.  Distance  between  any  point  on  the  ring  axis  and  the  center 
of  the  node  circle  ( L ). 

= Component  of  ^ for  n'th  buckling  mode  shape  ( L ). 
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R,  R'  ■*  Radii  of  curvature  of  the  ring  in  the  unstrained  and 
strained  states  respectively  ( L ). 

B * Curved  length  factor,  -»(ds/d«^)«  (df/<j<p) x]  (L). 

I am  Moment  of  inertia  of  the  section  of  the  ring  having  unit 
width  { iA  ). 


ds,  ds’  = Curved  length  of  an  element  in  the  undeformed  and 
deformed  states  respectively  ( L ). 

A ds  = Total  extension  of  the  ring  axis  ( L ). 

, yj  =*  Deflections  of  the  ring  axis  in  X-  and  Y-  directions 
respectively  ( L ). 

,db/  a Changes  of  length  of  dx  and  dy  respectively  ( L ). 
p =»  Uniform  external  pressure  per  unit  length  of  the  ring  of 
the  shell  ( MT~*  ). 


pQr  =s  Critical  pressure  of  a circular  ring  subjected  to  uniform 
external  pressure,  » (nl  - UEI/R*  , yj  « 2,  4, -*  • C MT“2  ). 
H,  V » Horizontal  and  vertical  components  of  force  per  unit  length 
( MT  -1  L ) . 


N -sw  Normal  foroe,  per  unit  length,  acting  through  the  centroidal 
axis  ( MT~2  L). 

Q «=.  Shearing  force,  per  unit  length,  parallel  to  the  cross 
section  ( MT  ~z  L ) . 

M =»  Bending  moment,  per  unit  length,  acting  at  any  point 
( MT"2  L2  ). 

M ^ =a  Component  of  M for  n’th  buckling  mode  shape  ( MT“2  L2  ). 
m*  =3  Static  moment,  per  unit  length,  due  to  external  loads 
only  { MT  ~2  L2  ) • 
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Ai&5  = Isorement  of  m5  resulting  from  the  deflection  of  the 
ring  axis  ( MT'2L2  J. 

(TMj<Tn—  Unit  fiber  stresses,  per  unit  length,  due  to  bending 
moment  and  normal  force  respectively  ( MT  L ) • 

(Add  )Jt  (Ad©  )x  — Rotations  of  cross  section:  due  to  the  exten- 
sion of  the  ring  axis  and  due  to  bending 
moment  respectively  ( Rad.). 

Add  =•  Total  rotation  of  cross  section  of  an  element  ( Rad.). 

A6  = Angle-change  from  a reference  point  to  any  point  ( Rad.  ). 
8 ~ Angle  between  the  tangent  and  X-axis  at  any  point  ( Rad.  ). 

=•  Angular  amplitude  of  a local  dent  in  a ring  ( Rad.  ). 

§=  Relati\e  amplitude  of  initial  irregularity  of  either  a ring 
or  a shell  ( Dimensionless). 

£ ss.  Unit  deformation  of  the  central  line  ( Dimensionless  )• 

e =*•  Elliptic ity  of  an  ellipse,  equal  to  1 minus  the  ratio  of 

minor  axis  to  major  axis  ( Dimensionless  ). 

n ss  Number  of  buckling  mode  shape  ( Dimensionless  ). 

^ =s  A constant,  =s  pRQ/EI  (Dimensionless  ). 

V = Poisson's  ratio  ( Dimensionless  ). 

F^"*  =*  Magnification  factor  for  n'th  buckling  mode  shape, 

=»  1/(1  - p/p  ) ( Dimensionless  ). 
cr . 

a,  , b„  ss  Coefficient  in  Fourier  Series  ( Dimensionless  )• 

D =*  Flexural  rigidity  of  the  seotion,  per  unit  length, of  a 
shell,  = Eh*/12(l  - V1  ) ( MT ” 2 L * ). 

K » A constant,  * (1  - V*  )/Eh  ( M 1 T2  ). 

A Curved  length  factor  in  longitudinal  direction,  » 

£l  + (d  jP/dx)2]  ( Dimensionless  ). 
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F =•  A constant,  =3  12(Ho/h  ) (Dimensionless  ) . 

x,<f  sa  Cylindrical  coordinates  ( L,  Rad.  ). 
m =.  A constant,  =•  TT  /L  ( L"1  ). 

2L  — Wave-length  of  the  dent  ( initial  imperfection  ) along 
the  longitudinal  axis  of  a shell  ( L ). 

X' , Y* , Z*  = Pressure  acting  on  the  shell  element  in  the 
directions  of  coordinate  axes  X,  Y,  and  Z 
respectively  ( MT“*  ). 

— Coefficients  of  deflection  configurations  w,  v, 
and  u of  the  middle  surface  of  the  shell 
( Dimensionless  ). 

Strains,  Displacements,  Stresses,  forces  and  moments  in  cylin- 
drical coordinates  folio*  the  notations  used  in  "Theory  of 
Plates  and  Shells"  by  S.  Timoshenko  (20). 
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